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ABSTRACT 

We present the Survey for High-z Absorption Red and Dead Sources (SHARDS), an ESO/GTC 
Large Program carried out with the OSIRIS instrument on the 10.4m Gran Telescopio Canarias 
(GTC). SHARDS is an ultra-deep optical spectro-photometric survey of the GOODS-N field covering 
130 arcmin 2 at wavelengths between 500 and 950 nm with 24 contiguous medium-band filters (pro- 
viding a spectral resolution R~50). The data reach an AB magnitude of 26.5 (at least at a 3cr level) 
with sub-arcsec seeing in all bands. SHARDS main goal is obtaining accurate physical properties of 
intermediate and high-z galaxies using well-sampled optical SEDs with sufficient spectral resolution 
to measure absorption and emission features, whose analysis will provide reliable stellar population 
and AGN parameters. Among the different populations of high-z galaxies, SHARDS principal targets 
are massive quiescent galaxies at z>l, whose existence is one of the major challenges of current hier- 
archical models of galaxy formation. In this paper, we outline the observational strategy and include 
a detailed discussion of the special reduction and calibration procedures which should be applied to 
the GTC/OSIRIS data. An assessment of the SHARDS data quality is also performed. We present 
science demonstration results about the detection and study of emission-line galaxies (star-forming 
objects and AGN) at z=0-5. We also analyze the SEDs for a sample of 27 quiescent massive galaxies 
with spectroscopic redshifts in the range 1.0<z<1.4. We discuss on the improvements introduced by 
the SHARDS dataset in the analysis of their star formation history and stellar properties. We discuss 
the systematics arising from the use of different stellar population libraries, typical in this kind of 
studies. Averaging the results from the different libraries, we find that the UV-to-MIR SEDs of the 
massive quiescent galaxies at z=1.0-1.4 are well described by an exponential decaying star formation 
history with scale t=100-200 Myr, age around 1.5-2.0 Gyr, solar or slightly sub-solar metallicity, and 
moderate extinction, A(V)^0.5 mag. We also find that galaxies with masses above M* are typically 
older than lighter galaxies, as expected in a downsizing scenario of galaxy formation. This trend 
is, however, model dependent, i.e., it is significantly more evident in the results obtained with some 
stellar population synthesis libraries and almost absent in others. 

Subject headings: galaxies: starburst — galaxies: photometry — galaxies: high-redshift — infrared: 
galaxies. 
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1. INTRODUCTION 

The current paradigm of galaxy formation estab- 
lishes that the baryons closely follow the evolution 
of the Cold Dark Matter (CDM) halos, which clus- 
ter and grow hierarchically as shown in cosmologi- 
cal simul ations and semi-ana lytical mod els (such as 
those in iSpringel et al.l 120051: see also iBaugh et al.1 

ISnmervillp kr PrimnrU 11 QQQl- iPnlp of nl I 190001- 



1998; ISomerville fc Primacfl 119991: iCole et all 
Somerville et al.ll2008t iRicciardelli fc Francesch ini 2010). 
In this scenario, star formation started within the cooling 
gas clouds in merging dark matter halos after a relatively 
slow early collapse regulated by feedback processes. This 
early star formation produced relatively small disk sys- 
tems that later merged and generated larger (i.e., more 
massi ve) spheroidal sy stems (see, e.g.. IKauffman n et al.1 
Il99l lEllis era!1l2000[ l200lh 

The global picture about the co-evolution of matter in 
the Universe (including all graviting components: CDM 
and baryons) is self-consistent and has been success- 
ful in reproducing and even predicting many observ- 
ables about galaxy evolution, especially at low redshift. 
Among the most relevant successes, we find the good 
comparison of models with the observed p ower spectrum 
of the Cosmic Microwave Background (|Spergel et al.l 
120071: iKomatsu et al.l 120111) or the Large Scale Struc- 
ture of the Universe (jPercival et al.l 12001b iColless et al.1 
120011 ). Also very convincing is the link between ob- 
servations and theoretical expectations such as the ex- 
istence and properties of the acoustic baryonic oscil- 
lations (lEisenstein et al.l 120051 : iGaztanaga et al.1 120091 : 
iPercival et al] 120101) . In addition, the hierarchical sce- 
nario for galaxy formation is also supported by the ob- 
servations o f galaxy mergers at different cosmological dis- 
tances (e.g., van Dokkum et aDfl999t Ivan DokkumH "2005: 



IMilev et all 12006). and the increase of the fraction of 
galaxies undergoing mergers as we move to higher red- 
shifts (among others. lLacev & Cole 1993; Lc Fevrc et al 
2000; IConselice et all 12003] iBell et al.1 120061: lLotz et a! 
2008', lLopez-Saniuan et al.N2009l ). 



However, the hierarchical picture contrasts with sev- 
eral observational evidences, especially at high redshift 
(z>l-2), where a more classical monolithic collapse is 
favored. This formation path was proposed 50 years ago 
to explain the origin of bulges such as the Milky Way's 
and spheroidal galaxies. This would be done through a 
free-fall rapid collapse causing the formation of the bulk 
of the stars in these systems in a short period of time. 
Later, the star formation is shut off by some quench- 
ing ph enomena, and the galaxy hence forth evolves pas- 
sively (jEggen et al.l 1962; La rson! 119741 ). This theory was 
largely abandoned due, first, to the compilation of ev- 
idence suppor ting that spheroidal gala xies suffer merg- 
ing episodes ()Toomre &: Toonirel 119721 ). Furthermore, 
globular clusters and the general stellar population in 
the MW present a relatively wide range of ages (e.g., 
iSearle fc Zinnl[l978l ). directly pointing out to the hierar- 
chical scenario. Eventually, the hierarchical picture was 
adopted instead of the monolithic collapse due to the 
high degree of success of the ACDM models and semi- 
analytic models mentioned above. 

Nevertheless, rapid early episodes of intense star for- 
mation are indeed consistent (although not uniquely) 
with observational facts in nearby galaxies, such as the 



dominant old stellar populations in bulges and ellip- 
ticals, their metallicity and a-elements enhancement, 
and the dynam ics an d shape of these s ystems (e.g., 
Vazdekis et all 119971 : IFaber et al.l 119971 : iTrager et all 
2000bl lal)" In addition, hierarchical models still present 
severe drawbacks in several aspects. The most chal- 
lenging observational facts for hierarchical models re- 
fer to the lightest and heaviest galaxies. Indeed, hi- 
erarchical models typically present a "missing satellite 
problem", i.e., they predict many mor e low-mass galax- 
ies than wh a t is actually observe d (see IKa uffman n et al] 
19931 IT999I: iKlypin et all Il99l iKravtsov et al.l 12004" 



Li & White 2009:~Quilis fc Truiilloll2012l) . On the bright 



end, models tend also to overpredict the number of mas 
sive galaxies observed in the local Universe, although 
they are getting closer to the observatio ns after tak- 
ing i n to account quenching mechanisms dCroton et al 
2001 iDe Lucia fc BlaizoH [20071: ISomerville et al.l 1200 ' 



Ricciardelli fc F ranceschini 20101 IGuo et al.ll2010D . 



The discrepancies between the predictions of current 
galaxy formation models based on the ACDM paradigm 
and the data are more obvious as we move to higher red- 
shifts. In the last 15 years, a wide variety of papers us- 
ing very heterogeneous data and methods have presented 
compelling evidence that the formati on of galaxies fol- 
lows a so-called downsizing scenario (ICowie et al.1 1 1996: 
i Heavens et al]|2004l: iGrazebrook et al 1120041: iBauer et al] 
20051: iPerez-Gonzalez et al l 120051: lArnouts et al.l 120071 : 
IPerez-Gonzalez et al.ll2008| ). In this theory, the most 
massive galaxies formed first in the history of the Uni- 
verse, thus having the oldest stellar populations seen to- 
day. The formation of less massive systems continued at 
lower redshifts. Downsizing implies that the bulk of the 
star formation in the most massive galaxies happened 
quick and stopped for some reason in early times. This 
also means that there should be massive passively evolv- 
ing galaxies at high redshift. This kind of objects have 
indeed been detecte d at redshifts around z~l-3 with a 
variety of techniqu e s (TYan et alj|2000atlFranx et al]|2003l: 
i Daddi et all 120041 : iPapovich et al.l 120061: ICimatti et"aLl 
2008). 

The finding of massive galaxies at z=l~3, some of 
them already evolving passively, is indeed extremely 
challenging for current models of galaxy formation 
based on the ACDM paradigm. Indeed, models pre- 
dict many less massive systems at high-z than observed 
(see, e.g.. IConselice et al.l 12007b iMarchesini et al.l 120091 : 
iHenriques et al.l 12012b iLeauthaud et al.l 120121 ). ~ In con- 
trast, the downsizing scenario contradicts, at least at first 
sight, the predictions of a hierarchical assembly of the 
stellar mass in galaxies, i.e., the most massive galaxies 
do not seem to be the result of multiple merg ers occurring 
in a extended period a l ong the Hubble tim e (jBaugh et al.l 
[1991 ICole et^l2000l IFaber et al.ll2007T ) . Still, a hierar- 
chical assembly with (maybe multiple) mergers occurring 
at high redshift between gas-rich systems (a process close 
in nature to a monolithic collapse) would be consistent 
with both the evidences for downsizing and the proper- 
tics of the dominant stellar populations seen in nearby 
spheroidal systems (e.g. lDekel et al.ll2009l ). 

From the observational point of view, our understand- 
ing of the processes involved in the early (z>l) assem- 
bly of galaxies (and also the evolution from the early 
Universe to the present) is still hampered by the sig- 



nificant (often systematic) uncertainties in our estima- 
tions of their physical properties. Our global picture 
of galaxy formation will only improve if we are able to 
get more robust estimations of some key properties of 
galaxies, such as the stellar masses, SFRs, and extinc- 
tions. Jointly with those, we of course need better es- 
timations of the distances to galaxies based on spectro- 
scopic or photometric redshifts, which can be used to 
relate the mentioned galaxy properties with other rel- 
evant parameters such as the environment. The im- 
provements in the determination of stellar masses and 
SFRs/extinctions will also mean a better estimation 



mation Histories, SFH (see, e. , 



Fontana et £01120061; iKriek et all 

Pforr et al.ll2012t lPacifici et al.lfeoTl i. Jointly with this 
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observational effort, models should also be improved, in- 
cluding better physics. For example, models are still to 
provide more certain emissivities of the stellar popula- 
tions in the rest-frame NIR, now affected by strong un- 
certainties due to limitations of about knowledge about 
the properties and importance of stellar evolutionary 
pha ses such as the thermally- pulsatin g TP-AGB phase 
(see lMarasto"i] [20051 : iKriek et al.l[2010h . The task of ob- 
taining more robust physical parameters of galaxies at 
cosmological distances is even more interesting for those 
massive galaxies which have already reached a quiescent 
state and are evolving passively at high-z, whose number 
densities and properties are the most demanding chal- 
lenges for current galaxy evolution models. The cosmo- 
logical importance of these systems is very high, since 
they most probably represent the early formation phases 
of present-day early-type galaxies. 

In this paper, we present the basics of the Survey for 
High-z Absorption Red and Dead Sources (SHARDS), 
an ESO/GTC Large Program awarded 180 hours of 
GTC/OSIRIS time during 2010-2013. This project 
consists of an ultra-deep (m<26.5 AB mag) imaging 
survey in 24 medium-band filters covering the wave- 
length range between 500 and 950 nm and targeting 
the GOODS-N field. The observations carried out 
by SHARDS allow to accurately determine the main 
properties of the stellar populations present in these 
galaxies through spectro-photometric data with a res- 
olution R~50, sufficient to measure absorption indices 
such as the D(4000) (e . g.. iBruzuall 119831: iBaloeh et al.l 
119991 iKauffmann et al.l l2003at IKriek et al.l 1201 lh or 
Mgr/y index (iSpinrad et al.lll997tlMcCarthv et al.ll200i 
Saracco et ail 120051: iDaddi et all l2005bl: iCimatti et al.l 



20081) . The analysis of these spectral features is a pow- 



erful method to constrain the solutions of stellar popula- 
tion synthesis models and to improve our estimations of 
parameters such as the age, SFH, mass, and extinction 
of galaxies at cosmological distances. 

SHARDS inherits the observational strategy of 
ast and on - going optica l surveys such as COMB017 
Wolf et all [20011 120031), t he COSMOS medium- 
band survey (illbert et al.l 120091 ) . ALHAMBRA 



dMoles et all 12008ft. a nd P AU/J-PAS (jBemtez et al.l 
I2009al lAbramo et al.l 120111 ). These projects have 
demonstrated the impact of large photometric datasets 
on our understandin g of the formation of galaxies (see , 
amon g many pap ers, iBell et all 120041: iFaber et al.l 120071: 
iWolf et all 12001 iRix et al.l l2004t feorch et al.l boOS 



Caputi et all 120061: iScoville et all l2007t iBem'tez et al.l 
2009b nCardamone et al.l 120101: iWhitaker et al. I 12011ft . 
SHARDS intends to be a step forward from these 
surveys in terms of depth, spectral resolution, and data 
quality. Our survey prioritizes the detailed study of 
the faintest galaxies at the highest redshifts over the 
analysis of closer galaxy populations and the Large Scale 
Structure at intermediate redshift, and thus focuses 
on a smaller area than the surveys mentioned above. 
Indeed, SHARDS was planned to reach up to 3 mag 
fainter than those surveys, uses typically twice the 
number of filters in the same wavelength range (i.e., 
our spectral resolution is better), and was obtained in 
excellent (sub-arcsec) seeing conditions with a 10m class 
telescope. In contrast, it covers a fraction of the area 
surveyed by other projects. 

In this paper, we present the main technical charac- 
teristics of the survey in Section O including a thorough 
discussion of the reduction and calibration procedures 
in Section [3j Next, we present our science verification 
results about emission-line and absorption systems. In 
Sections 14.11 and 14 . 2 [ we discuss about the ability of the 
SHARDS data to select and study emission-line sources 
(star- forming galaxies and AGN) at intermediate (z<l) 
and high redshifts (up to z^5 and beyond). In Sec- 
tion [SI we present detailed spectral energy distributions 
of massive quiescent galaxies at z>l, and demonstrate 
the power of our spectro-photometric data to analyze the 
stellar populations in this kind of object through a de- 
tailed comparison with stellar population synthesis mod- 
els. 

Throughout this paper we use AB magnitudes. We 
adopt the cosmology H = 70 km s _1 Mpc _1 , Q m = 0.3, 
and f2 A = 0.7. 

2. SURVEY DESCRIPTION 

The Survey for High-z Absorption Red and Dead 
Sources is a medium-band optical survey currently be- 
ing carried out with the Spanish 10.4m telescope, Gran 
Telescopio Canarias (GTC), and its OSIRIS instrument. 
SHARDS was approved in 2010 as an ESO/GTC Large 
Program and awarded with 180 hours of observing time 
to obtain data through 24 contiguous medium-band fil- 
ters covering the wavelength range between 500 and 950 
nm. The survey targets the GOODS-N field, covering 
most of the area observed with the HST/ACS instru- 
ment. SHARDS was conceived to study in detail the 
properties of the stellar populations in 0<z<4 galax- 
ies (and beyond), focusing in the major goal of analyz- 
ing quiescent massive galaxies at z= 1.0-2. 5. To achieve 
this goal, the survey was planned to obtain photomet- 
ric data in as many filters as necessary to cover the 
entire optical window with enough spectral resolution 
to be able to measure reliably absorption indices which 
could be used to perform a detailed and robust stellar 
population synthesis. Indices such as D(4000) or Mgyy 
use spectral windows 10-20 nm widths. To obtain this 
spectral resolution (or better) at z>l, we imposed fil- 
ter widths of approximately 15-20 nm. A compromise 
between spectral resolution, depth, and manufacturing 
limitations was adopted and the survey was planned to 
use 17 nm wide filters. The bright night sky at wave- 
lengths beyond ^800 nm imposes that filters as narrow 
as 17 nm would not reach the needed magnitude limit to 
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Fig. 1. — Observing strategy of SHARDS. The figure shows a 
scheme of the transmission curves for the 24 SHARDS medium- 
band filters (black and red lines at the bottom, with names giving 
nominal central wavelengths and widths, as written on the leg- 
end). This filter set was designed to probe the optical wavelength 
range contiguously between 500 and 950 nm with filters of width 
FWHM~17 nm (spectral resolution R~50). The observing strat- 
egy of SHARDS was devised to identify quiescent galaxies at z>l. 
We show how the stacked spectrum of the 13 spectroscopically 
confirmed quiescent massive galaxies at 1.4<z<2.0 in GOODS-S 
HCimatti et al.ll20Qgl. the stack adds up 480 hours of VLT time) 
would look like at 4 different redshifts. The sky spectrum is also 
depicted in green. Using 24 filters, SHARDS probes the prominent 
absorption feature placed at A=265-295 nm, distinctive of passively 
evolving galaxies (marked in red) with a resolution R~50, and it is 
able to measure its strength with a spectro-photometric technique 
to accurately determine stellar age and distances for individual 
high-z galaxies. The 900 nm atmospheric window is not covered by 
our survey, since it will be cover ed by another approved GTC Guar- 
antccd Time project (OTELO, Ccpa ct al. 2011; Lara-Lopez ct al. 

Eon]). 



study high-z sources, so our 2 reddest filters are twice as 
wide. SHARDS was conceived to reach the typical mag- 
nitudes of sub-L* galaxies at z>l for every single filter, 
so the goal was to obtain a depth of 26.5-27.0 mag at the 
3rj level with sub-arcsec seeing. The final manufactur- 
ing process resulted in an average width for our filter set 
around 15-16 nm, except for the 2 reddest filters, which 
have a width of 33-35 nm. Figure [1] shows the layout of 
the filter set and the observational strategy of SHARDS. 
The characteristics of the filters are given in Table [IJ 
jointly with other details of the SHARDS data. 

SHARDS is being carried out in the GOODS-N field, 
one of the most targeted areas of the sky at all wave- 
lengths. Virtually all the deep region covered by ACS 
is being surveyed by SHARDS using two OSIRIS point- 
ings, summing up a total surveyed area of ^130 arcmin 2 
(see Figure [2]). The central J2000 coordinates of 
the two pointings are a=12:37:18.9 <5=+62:17:03 and 
a=12:36:33.3 <5=+62:ll:39. A position angle of 45° was 
used for our imaging data in order to cover the GOODS- 
N region more efficiently. 

The multi-wavelength dataset available in GOODS- 
N is extensive, ranging from an ultra-deep X-ray ex- 
posure to the deepest da ta in the MIR/FI R with 
surveys such as FIDEL dFraver et all I2006D. PEP 
(|Lutz et al.1 EoH , HerMES (IQliyer et all 120101). or 
Herschel-GOODS (lElbaz et al.ll2011l) as wel l as multiple 
spectroscopic redshifts (jWirth e t al. 2004; Cowi e et all 
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Fig. 2. — Footprint of the SHARDS data overimposed on the 
ACS images for the GOODS-N field. The footprints of the HST 
grism surveys carried out with ACS (PEARS) and WFC3 (PI: 
Weiner) and the CANDELS coverage of GOODS-N are also shown. 
SHARDS covers a total surveyed area of ~130 arcmin 2 divided in 
two pointings. 
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[200l iReddv et al l [2001. a few of them at z>1.5 
(|Reddv et al.1 120051 : iBarger et al.ll2008| ). Very comple- 



mentary to the SHARDS data, GOODS-N has been ob- 
served by HST with ACS and WFC3 (see Figure 
providing slitless, intermediate resolution spectroscopy 
in the optical (thr ough the G800L grism; PEARS , 
IKummel etail 120091: see also IPirzkal et aLI 12004 I2009j) 
and NIR (G141; PI: B. Weiner). In addition, the avail- 
ability of the sky deepest IRAC ([3.6]<26.0 mag) and 
MIPS [F5 CT (24)>30 /LtJy] observations ensures the detec- 
tion of the rest-frame NIR/MIR emission of the galax- 
ies and allow us to robustly estimate stel l ar masses and 
SFRs (jPerez-Gonzalez et all I2005L 120081 : iCimatti et all 
120081 ). This wealth of data converts GOODS-N in one 
of the two best fields (the other being GOODS-S) for the 
study of the first galaxies and their evolution, and the 
best in the northern sky hemisphere. 

At the time of publication of this paper, 75% of the 
SHARDS data have already been taken. Table Q] presents 
the data characteristics for the observed filters in the two 
SHARDS pointings. 

3. DATA REDUCTION AND CALIBRATION 

3.1. Reduction pipeline 

The GTC/OSIRIS instrument (IC^IMol) is an op- 
tical imager equipped with 2 CCDs, which cover a to- 
tal (usable) field-of-view of 7.4' x 8.4'. There is a gap 
between the 2 detectors, with size 10-12" as measured 
directly in our images in different locations along the 
gap and after masking non-useful pixels and correcting 
for distortions. OSIRIS presents special characteristics 
which should be taken into account in the data reduc- 
tion procedure. Consequently, apart from the standard 
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TABLE 1 

Characteristics of the SHARDS filter set and observations (before 2012A). 



Filter 


CWL 


Width 


A 


B 


X 


Y 


rms 


AZP 


m3 


(7 


m 75% 


seeing 


















PI 


P2 


PI 


P2 


PI 


P2 


PI 


P2 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(9) 


(10) 


(10) 


(11) 


(11) 


(12) 


(12) 


F619W17 


618.9 


15.1 


623.14 


-2.402e-6 


-202 


985 


0.10 


0.077 




26.89 




27.22 




0.85 




F636W17 


638.4 


15.4 


641.37 


-2.587e-6 


-116 


986 


0.10 


0.066 


0.072 


26.78 


26.70 


27.15 


27.17 


0.79 


0.92 


F653W17 


653.1 


14.8 


656.01 


-2.634e-6 


-151 


999 


0.10 


0.063 


0.065 


26.91 


27.07 


27.17 


27.15 


0.98 


1.00 


F670W17 


668.4 


15.3 


671.86 


-2.600e-6 


-183 


1037 


0.10 


0.062 


0.064 


26.64 


26.76 


27.12 


27.18 


0.79 


1.07 


F687W17 


688.2 


15.3 


690.50 


-2.674e-6 


-186 


983 


0.10 


0.050 


0.055 


27.04 


26.85 


27.10 


27.09 


0.84 


0.93 


F704W17 


704.5 


17.1 


707.78 


-2.723e-6 


-209 


1028 


0.10 


0.055 


0.063 


26.71 


26.63 


27.01 


26.95 


0.89 


0.92 


F721W17 


720.2 


18.2 


723.12 


-2.969e-6 


-94 


960 


0.10 


0.056 


0.060 


26.60 


26.51 


26.97 


26.96 


0.93 


1.02 


F738W17 


737.8 


15.0 


741.80 


-2.411e-6 


-328 


1050 


0.10 


0.050 


0.061 


26.45 


26.25 


26.95 


26.87 


0.86 


0.90 


F755W17 


754.5 


14.8 


758.12 


-2.660e-6 


-228 


1034 


0.10 


0.050 


0.055 


26.69 


26.37 


26.93 


26.91 


0.92 


0.93 


F772W17 


770.9 


15.4 


774.62 


-2.929c-6 


-122 


1026 


0.10 


0.054 


0.053 


26.54 


26.34 


26.90 


26.83 


0.94 


1.04 


F789W17 


789.0 


15.5 


791.22 


-3.087e-6 


-123 


994 


0.10 


0.054 


0.056 


26.22 


26.02 


26.84 


26.63 


0.97 


0.92 


F806W17 


805.6 


15.6 


809.42 


-2.939c-6 


-200 


933 


0.10 


0.051 


0.050 


26.39 


26.38 


26.82 


26.77 


0.96 


0.99 


F823W17 


825.4 


14.7 


829.15 


-3.055e-6 


-153 


888 


0.10 


0.049 


0.049 


26.59 


26.65 


26.91 


26.89 


0.82 


0.91 


F840W17 


840.0 


15.4 


843.51 


-3.103e-6 


-237 


992 


0.10 


0.057 


0.056 


26.13 


26.19 


26.79 


26.74 


0.88 


0.91 


F857W17 


856.4 


15.8 


859.97 


-2.892e-6 


-249 


1002 


0.10 


0.050 


0.051 


26.84 


26.23 


26.88 


26.63 


0.73 


0.95 


F883W35 


880.3 


31.7 


885.33 


-2.889e-6 


-285 


978 


0.10 


0.065 


0.057 


26.06 


26.06 


26.64 


26.58 


0.93 


1.02 



Note. — (1) Filter name. (2) Central wavelength (in nm) of the filter for angle of incidence AOI— 10.5° (approximately that for the center of the 
FOV). (3) Filter width (in nm). (4) Coefficient A (in nm) for CWL calibration (from equation^. (5) Coefficient B (in pixel" 2 ) for CWL calibration 
(from equation 0. (6) X coordinate for the OSIRIS optical center (in pixels). (7) Y coordinate for the OSIRIS optical center (in pixels). (8) RMS 
of the CWL calibration (in nm). (9) Zcropoint uncertainty for pointings 1 and 2. (10) Sensitivity limit at 3-er level (AB mag) for pointings 1 and 2. 
(11) Third quartil of the magnitude distribution (AB mag) for pointings 1 and 2. (12) Average seeing (in arcsec) for pointings 1 and 2. 



reduction steps (bias subtraction and flat-fielding), our 
custom OSIRIS pipeline includes the following additional 
steps: illumination correction, background gradient sub- 
traction, fringing removal, WCS alignment taking into 
account field distortions, and 2-dimensional calibration 
of the photometric pass-band and zeropoint. In what 
follows, we elaborate on these steps and describe in de- 
tail the non-standard parts of the reduction. 

Bias frames were taken each night. These calibration 
data were found to be very stable and uniform, with typi- 
cal variations of less than 3% from pixel to pixel and from 
night to night. Dark current was found to be negligible 
in all our data. 

Sky flats for each filter in the SHARDS dataset were 
taken at sunset and dawn every night. These flats pre- 
sented strong gradients across the field (typically, 10-50% 
differences in brightness from one edge of the FOV to an- 
other), with significant spatial and temporal brightness 
variations. These gradients were also seen in science data 
during both dark and (especially) gray night time. 

After analyzing night and day-time data, we concluded 
that these spatial variations were mainly due to the 
special characteristics of the OSIRIS instrument, which 
operates off-axis for our observational setup (medium- 
and broad-band imaging). Indeed, light rays reach the 
OSIRIS detector (when not using the tunable filters) in a 
wide range of angles of incidence (AOI): -2°<AOI<22°. 
Given the typical dependence with the light ray's AOI 
of the central wavelength (CWL) of the passband for in- 
terference filters, and the medium- width of our filter set, 
our flat-field frames were subject to the spectral features 
of the sky spectrum, and their temporal variations. 

The optical axis is outside the OSIRIS FOV when using 
a single broad- or medium-band filter. For this reason, 
a gradient is observed in our data (flat-field and science 
images) with a symmetry around an horizontal line ap- 
proximately dissecting in equal parts the FOV in the ver- 



tical direction. The gradient follows a radial profile cen- 
tered in the optical axis. It presents varying structures 
of different brightnesses whose position and strength de- 
pend on time and on the filter. The brightest structures 
are located where strong sky emission lines (or bands) 
go through the pass-band for each physical filter as its 
CWL is changing along the FOV. 

This effect directly related to the special characteris- 
tics of OSIRIS and our instrumental setup means that it 
is very difficult to find a spatially constant light source to 
take flat- fields, i.e., a significant part of the structure of 
the flat-field frames is linked to the sky spectrum. Thus, 
our pipeline included an illumination correction to get rid 
of this effect as much as possible. This illumination cor- 
rection was carried out by comparing the flat-field images 
taken for our medium-band filters with those acquired 
through broad-band pass-bands (typically, r- and z-band 
filters). For these broad-band filters, the sky spectrum 
is averaged in a spectral range which is wide enough to 
prevent strong spatial variations, and are also very con- 
stant in time. Indeed, the super-skyflat provided by the 
observatory and built with thousands of frames shows 
variations of less than 5% along the FOV and from night 
to night. 

We used the ratio of our medium-band flat-fields to 
the broad-band flat-fields to correct the SHARDS data 
for the illumination effect. These ratios were smoothed 
using a 3 th -order spline. 

After applying the flat-field and illumination correc- 
tions to the science data, these images presented a highly 
symmetric sky gradient, which was subtracted with a 
median filtering and spline interpolation, after masking 
objects. In order to avoid the effect of the wings of the 
objects in this sky determination, we increased the exten- 
sion of the sourc es (typically by a factor of 2 in Kron ra- 
dius, [Kron 1980) taking into account their brightnesses, 
and we also masked out faint objects detected in a pre- 



liminary mosaic constructed for each filter adding all the 
available observations. 

The data for the filters whose CWLs are redder than 
^700 nm presented some fringing (typically with an in- 
tensity below 1% of the background). We removed this 
additive effect using the RMFRINGE task in IRAE0. 

Before stacking all the data together, we calibrated the 
World Coordinate System (WCS) for each image using 
the positions of several hundred objects for each CCD, 
which were cross-correlated with a sample of galaxies de- 
tected in the Subaru i?-band image of G OODS-N. Fo r 
this task, we used the SCamp program ([Bertinl [2006). 
and the WCS utilities in IRAF, obtaining for each frame 
an undistorted image remapped to a TAN-SIP coordi- 
nate system. The distortion in individual frames reaches 
several arcseconds in the edges, and the final WCS cal- 
ibrated images present a typical position uncertainty 
around 0.1" throughout the FOV. 

Finally, the SHARDS pipeline stacks together all the 
data for a given filter with a sigma-clipping algorithm 
to get rid of cosmic rays and artifacts, producing final 
mosaics and exposure maps with the same WCS for all 
SHARDS filters and for each of our 2 pointings covering 
the GOODS-N field. 

3.2. Calibration procedures 

Given the special characteristics of the OSIRIS in- 
strument at GTC, within each single frame taken with 
a given physical filter, each pixel sees a different pass- 
band. This has a strong effect on the images taken with 
medium-band filters, where the shift of the central wave- 
length of the actual pass-band seen by different parts 
of the detector produces sky gradients and significant 
difficulties in obtaining flat-field calibration images, as 
explained in the previous section. To overcome this is- 
sue, we performed a detailed calibration of the effective 
pass-band and photometric zero-point as a function of 
the position in the FOV. Although part of the spatial 
variation of the photometric zeropoint can be removed 
by applying an illumination correction built with broad- 
band flat-field images, a small effect can remain in the 
images. To account for this, during our calibration pro- 
cedure we considered that the zeropoint could vary along 
the image. The calibration method was devised, conse- 
quently, to determine the spatial variation of the zero- 
point. Note that after the calibration routine, we can 
construct a photometric catalog for each physical filter, 
but within this catalog, the flux measurements for each 
detected source refer to a different pass-band. This is 
quite different from a standard photometric catalog for 
a given physical filter, where all sources share the same 
pass-band. 

We describe the procedures to carry out the calibra- 
tion of the SHARDS data in the following subsections, 
starting from the CWL calibration and following with 
the absolute flux calibration. 

3.3. Pass-band central wavelength (CWL) calibration 

19 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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Fig. 3. — Top: Calibration of the transmission curve of the 
F687W17 filter as a function of angle-of-incidence, obtained with 
laboratory data. The central wavelengths of each pass-band are 
marked with vertical lines. Bottom: Calibration of the spatial vari- 
ation of the CWL of the F687W17 filter along the OSIRIS FOV. 
This variation is symmetric and only depends on the distance to 
the optical axis. The data are fitted with the function in Equa- 
tion [l] (green line), leaving the position of the optical axis as a 
fitting parameter. The best fit between CWL and distance to the 
optical axis in arcmin is given in the plot, jointly with its rms. 

The central wavelength of the pass-band seen by dif- 
ferent parts of the OSIRIS detector for a given physical 
filter varies along the FOV. This dependence was cal- 
ibrated with day-time imaging and spectroscopic data 
taken in laboratory and at GTC. In the lab, we cali- 
brated the transmission curve for each filter as a function 
of the angle-of-incidence of the incoming beam. At the 
telescope, once the given filter was mounted in OSIRIS, 
we took spectroscopic data through a special mask with 
105 pinholes homogeneously covering the entire FOV. 
The positions of those pinholes were determined with 
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images taken through the mask and using one of our fil- 
ters. Then, we also took spectra of these pinholes with 
the R1000B and R1000R grisms (which cover the entire 
wavelength range probed by SHARDS) through our fil- 
ters, calibrating them with Ne, HgAr, and Xe arcs. Typi- 
cal uncertainties in the wavelength solution were smaller 
than 0.01 nm. These data were used to measure the 
transmission curve at each pinhole position. We char- 
acterized each curve using two parameters: the CWL 
and the width of the pass-band. The shape and width 
of the pass-bands are relatively independent of the po- 
sition in the FOV (i.e., all pinholes show similar values 
within 0.2 nm). However, the CWL presents a signif- 
icant symmetric variation around the optical axis, lo- 
cated to the left of the FOV. The CWLs for all the pin- 
holes were fitted with a function depending on the square 
of the distance to the optical axis (r 2 ), leaving also the 
position of the optical ax is as a parameter to fit (see 
iMendez-Abreu et al.l 120111 for similar calibration proce- 
dure, but for the OSIRIS red tunable filter). An example 
of the calibration for one of the SHARDS filters is shown 
in Figure [3 The data were fitted with the following func- 
tion: 

CWL(X, Y) = A + Bx[(X- X ) 2 + (Y - Y ) 2 } (1) 

where X and Y are the positions (in pixels) in the OSIRIS 
FOV, and Xq and Yq are the position of the optical axis, 
combined all of them to give a distance to the optical 
axis in pixels. Xq, Yq, A, and B were measured for each 
one of the SHARDS filters by fitting the data described 
above and shown in Figure [3] for filter F687W17. When 
applied to the actual science data, the optical axis posi- 
tion was converted to RA and DEC, and the distances 
were measured in arcmin using the nominal pixel size of 
our data, 0.251 arcsec/pixel. 

The fitting coefficients in Equation [T] for each one of the 
SHARDS filters (with available data so far) are given in 
Table [TJ We want to remark that the four parameters 
in Equation [T] are highly correlated. Consequently, the 
differences seen from filter to filter may not be closely 
related to real differences in the position of the optical 
axis or CWL variation coefficients. In any case, we are 
only interested in recovering CWL values for any position 
within the FOV, and even assuming large (correlated) 
uncertainties for the four parameters, the results given 
in Tabled] do provide very accurate CWLs. In addition, 
we tested the repeatability of these measurements and 
no significant changes in this calibration were detected 
for observations taken in different nights. 

3.4. Photometric calibration 

The significant variation of the passband seen by each 
point of the detector as a function of the position in the 
FOV implies a complex behavior of the absolute photo- 
metric calibration of the SHARDS images. Moreover, the 
effects of this CWL variation on the construction of the 
flat-field (explained in Section 13. 1[) may also affect the 
flux calibration of the final mosaics. To cope with these 
issues, we developed a special flux calibration procedure, 
aimed at determining the zeropoint of the SHARDS mo- 
saics in each filter as a function of position in the im- 
age. Note that the behavior of the passband is symmet- 
ric around the optical axis and each position in the FOV 



is characterized by a central wavelength (as explained in 
Section 13.31 and defined by Equation [1]) . Keeping this in 
mind, we work with zeropoint variations as a function 
of CWL, instead of directly relating the calibration to a 
position in the images. 

The flux calibration of the SHARDS mosaics was per- 
formed by comparing the measured photometry in our 
images with spectroscopic data for several sources in the 
field (most of them being galaxies). Given the need for 
spectroscopic calibration data for sources covering the 
whole FOV, direct observations taken with GTC/OSIRIS 
were unaffordably expensive in terms of observing time, 
and even impossible since OSIRIS currently has limited 
spectroscopic capabilities (only long-slit). We therefore 
used spectroscopic data taken with other telescopes. In 
any case, a spectro-photometric standard star was ob- 
served with OSIRIS with a grism and each one of our 
filters, and we checked that our main calibration pro- 
cedure (described below) is consistent with these (very 
limited) observations. 

For our main calibration procedure, synthetic fluxes 
were obtained by convolving the spectra for sources dis- 
tributed all along the FOV with the appropriate trans- 
mission curve seen by that source, according to our cali- 
bration of the passband (in terms of the position in the 
FOV, see Section 13. 3p . The spectra used in this cali- 
bration method were taken from 2 different sources: (1) 
the HST ACS grism data from the PEARS project EJ 
and (2) the spectroscopic data released b y the Team 
Keck Treasury Re dshift Survey (TK RS IWirth et all 
l200l and DEEP3 (|Cooper et al.ll20Pl taken with the 
Keck/DEIMOS instrument. We note that the TKRS and 
DEEP3 spectra are not flux calibrated, they are not com- 
pletely flat- fielded in the spectral direction (see below), 
and are subject to the effect of strong sky emission lines 
and telluric absorption bands in certain spectral regions. 
For these reasons, our primary calibrator was the PEARS 
dataset of HST/ ACS spectra. The Keck spectra were 
used as a consistency check for the primary calibration. 
In addition, we also performed another test of the calibra- 
tion based on synthetic magnitudes obtained from stellar 
population models fitting the broad-band photometry for 
the galaxies with a published spectroscopic redshift and 
com piled in the Rainbow Cosmological Surveys dat abase 
(see lPerez-Gonzalez et al.ll20 08: Barr o et al.ll20 11 al ibi). A 
finer calibration will be carried out in the future based on 
high SNR (SNR>10) spectra for bright objects which are 
currently being obtained at TNG and the CAHA 3.5m 
telescope. The several methods to calibrate the data are 
explained in the following sub-sections. 

3.4.1. Primary calibrator: HST/ACS grism spectroscopic 

data 

The photometric calibration procedure starts with a 
crude determination of the zeropoint based on the com- 
parison of photometry in the SHARDS images with R- 
or /-band fluxes. This comparison provides a first es- 
timation of the conversion from counts to AB magni- 
tudes (one value for a given image). Then, we use the 
HST/ ACS spectra to improve this calibration and check 

20 Downloaded from the da tabase at 

http:/ /archive.stsci.edu/prepds/pears/ See also IFerreras et at] 
(2009). 
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Fig. 4. — Examples of the HST/ACS spectra (from the PEARS data described in Ferrcras ct al. 2009 and the PEARS database) used 
in the calibration of the SHARDS data. The ACS spectra are shown in black with the gray shaded area depicting the uncertainties. The 
spectra are scaled to the ACS viz photometry, shown in green, including also uncertainties and filter widths. The convolutions of the 
spectra with the viz filter transmission curves are shown with squares. The pass-band for the F687W17 filter seen by each galaxy is shown 
in blue, and the convolution with the HST spectra is marked with a red star, jointly with the convolutions for the other SHARDS filters 
observed so far (no transmission curves are depicted for the sake of clarity). We show the SHARDS and i-band ACS images of the galaxy 
(5" X 5") and the color b etween them as a function of aperture radius. In each plot, the best elliptical aperture determined by sextractor 
IIBertin & Arnouts 1993) f° r the ACS data is shown in green, the best aperture for the SHARDS image is plotted in red, and the circular 
aperture used to calibrate the SHARDS data is shown in black. The panel in the last row shows a galaxy whose photometry (and probably 
the spectra) are affected by contamination from nearby sources. This kind of source was excluded from our calibration procedure. They 
are marked with gray symbols in Figure 151 
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CWL [nm] 

Fig. 5.— Calibration of the SHARDS data for filter F636W17 
based on the HST/ACS spectra. The difference between the AB 
magnitude obtained from the spectra and the magnitude measured 
in the SHARDS images is plotted as a function of CWL of the 
pass-band seen by the gala xy (i .e., vs. position in the field-of- 
view, as described in Section 13.31 and defined by Equation [TJ. The 
points show the calibration obtained from individual galaxies for 
the best photometric aperture, i.e., that enclosing the entire object 
and minimizing the uncertainties, the effects of seeing, and the 
contamination from close neighbors. Error bars are estimated from 
the grism flux uncertainties and the rms of the comparison of ACS 
spectra and broad- band photometry in the bviz filters (0.4 mag). 
Objects whose photometry is affected by nearby sources or have 
low-SNR spectra are plotted with gray symbols. Color solid lines 
show the median behavior for the offset calculated with different 
photometric aperture radii. The final calibration is shown with 
a green dashed line, resulting from a polynomial fit to the best 
aperture data. 



the dependence on the position along the FOV. We start 
by scaling the ACS spectra to the broad-band photome- 
try in the viz bands from ACS. The rms of the compar- 
ison of broad-band photometry and the synthetic fluxes 
obtained with the grism spectra is 0.38 mag, i.e., this is 
the typical accuracy of the spectroscopic data for each 
galaxy. The flux measurement in the ACS bands is car- 
ried out within a large enough aperture so it can be com- 
pared with photometry in the SHARDS images with neg- 
ligible seeing effects and also avoiding (as much as pos- 
sible) the contamination from nearby sources. Typically, 
apertures with radii larger than 1" were used. The pro- 
cedure basics are outlined in Figure H) the spectra are 
scaled to the broad-band photometry, we measure the 
flux in the SHARDS pass-band corresponding to the po- 
sition of the galaxy in the FOV by convolving the spectra 
with the appropriate transmission curve, and we compare 
with the flux in counts measured in the SHARDS image, 
obtaining a conversion from counts to flux densities or 
AB magnitudes as a function of position in the FOV. 

Figure [3] shows the results of the calibration based on 
HST/ACS spectra. We plot the offset between the cal- 
ibration obtained from the analysis of the HST spectra 
and the preliminary calibration obtained with the com- 
parison with ii-band fluxes as a function of the CWL 



of the passband seen by each galaxy. In this example, 
we show the results for the observations of one of the 
SHARDS pointings through the F636W17 filter, where 
we have more than 1,500 ACS spectra to compare with. 
For the second pointing, PEARS only covered a part 
of it, and the number of available ACS spectra is ~900. 
The derived differential calibration as a function of wave- 
length (i.e., position) was applied to the images, so the 
final mosaics have a constant zero-point throughout the 
whole image. The typical scatter around the average cal- 
ibration was below 0.1 mag. Note that many spectra pro- 
vided very large zeropoint offsets. We visually inspected 
all spectra with offsets above 0.5 mag and virtually all of 
them presented nearby bright objects which were most 
probably contaminating the spectra (e.g., the last source 
in Figure 01). This effect was confirmed by comparing the 
scale factor of the spectra to the broad-band photometry. 
In spectra affected by contamination by nearby objects, 
the three bands used in the scaling provided very dif- 
ferent factors, with the scatter being considerably larger 
than the typical value (~0.4 mag). Note also that we did 
not establish any magnitude cut in the HST/ACS data, 
and a significant fraction of the spectra correspond to 
very faint objects and count with large uncertainties. 

3.4.2. Calibration test: ground-based spectroscopic data 

The calibration obtained with the HST data was 
checked with other spectra taken from the literature. 
We found nearly 1,000 spectra for sources detected by 
SHARDS in the TKRS database and the DEEP3 releas e 
in GOODS-N (IWirth et al.l 12001 ICooper et al1l20Tl . 
The spectra were all taken with the DEIMOS instru- 
ment on Keck with typical exposure times of 1-2 hours 
through R~300-600 grisms. Given that these surveys 
were mainly interested in spectroscopic redshifts based 
on emission-lines, the typical SNR in the continuum is 
low (SNR^l at the original spectral resolution). In ad- 
dition, the spectra are not flux calibrated and they typ- 
ically show a decrease of flux to the blue that seems 
to arise from the lack of an accurate spectral flat cor- 
rection (see Figure [6]). Moreover, these spectra were 
taken through a mask with slit widths of 1", so aper- 
ture effects may also be significant, i.e., the spectro- 
scopic data give information about the central part of 
the galaxies and a comparison with seeing limited pho- 
tometry can be biased. For these reasons, we just used 
these spectra as a test for our basic absolute flux cali- 
bration based on HST data. Prior to this test, we had 
to correct the individual ground-based spectra with an 
average spectral flat and flux calibrate them. Both pro- 
cedures were carried out using broad-band photometry 
and fits to this photometry with stellar population syn- 
thesis mo dels from the Rainbow Co smological Surveys 
database (Perez-Gonzalez et al. 2008). 

Figure [5] shows two typical spectra in GOODS-N ex- 
tracted from the ground-based spectroscopic sample. Af- 
ter binning the spectrum to a resolution of 2 nm, we 
carried out a preliminary flux calibration of the spec- 
trum based on the photometric broad-band data around 
800 nm (blue stars in the right panels of the figure) . An 
inspection of the resulting calibrated spectra (thick black 
lines in Figure H]) revealed that the spectroscopic data 
were not corrected with an adequate spectral flat, result- 
ing in a large and increasing flux difference between the 
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Fig. 6. — Two examples of the ground- based spectra used to test the absolute flux calibration of the SHARDS data. For each galaxy, 
the original spectrum (in arbitrary F u units and binned to achieve a resolution of 2 nm) is plotted on the left, indicati ng the name of the 
source in the Rain bow Cosmological Surveys database (Perez-Gonzalez et al. 2008) the TKRS and DEEP3 catalogs IWirth et al.ll2004l ; 
ICooper et al. 2011) and the spectroscopic redshift. Both on the left and right panels, we mark with a shaded area the regions where strong 
telluric absorption bands are located. On the right, we show with a black thick line the original spectrum after flux scaling it to match 
the broad-band photometry, depicted with blue stars (only the reddest points were used in this first scaling, given that they provided a 
better absolute calibration) and gray symbols (for the bluest bands which are more affected by the spectral flat and atmospheric extinction 
corrections). The black thin line shows the same spectrum after applying an average spectral flat obtained as explained in the text. 
The green line shows the final ground-based spectrum used for the calibration of SHARDS data, obtained after applying an atmospheric 
extinction correction calculated through a comparison with a stellar population synthesis model fitting the broad-band photometric data, 
shown in red in the plots. Black squares and circles, and green circles show the convolution of the previously described spectra (respectively: 
the original after flux calibration, the one after further calibration with a spectral flat-fielding, and the final one after carrying out also an 
extinction correction) with the transmission curves for the broad-band filters. 



spectra and the broad-band photometric data at shorter 
wavelengths (see comparison between black thick lines 
and the photometry in Figure [5]). To correct for this ef- 
fect, we built a master spectral fiat by averaging the ra- 
tio between the spectra and stellar population synthesis 
models fitting the broad-band photometry (red lines in 
Figure [6]) for all sources in the spectroscopic sample. Af- 
ter applying this correction, the comparison between the 
spectra (thin black lines in Figure[B]) and the photometry 
was better, but still most sources presented a dimmer flux 
level in the spectra when compared with the photometry 
for the bluest wavelengths. We identified this flux differ- 
ence as an effect of differential extinction, which was not 



fully taken into account in the average spectral flat. To 
account for this, we applied a final correction based on a 
typical extinction curve at Mauna Kea scaled to explain 
the difference between the spectra and the broad-band 
photometry (especially in the bluest bands). This final 
spectrum (green lines in Figure [5]) was then convolved 
with the SHARDS filter pass-band (the appropriate one 
according to the position of the galaxy in the FOV) , pro- 
viding another independent calibration for our data. 

3.4.3. Calibration test: synthetic photometry 

Finally, the fluxes obtained in the SHARDS bands were 
compared directly with synthetic magnitudes obtained 
by convolving the SHARDS filter pass-bands with stellar 
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population synthesis models fitting the broad-band pho- 
tometry for each source. Note that this calibration test 
is not completely independent of the one carried out with 
the ground-based spectra since those spectra were par- 
tially calibrated with the stellar population fits. More- 
over, the comparison with models for individual galax- 
ies may be biased due to the presence of emission-lines, 
which are not constrained by the broad-band photome- 
try. However, because of the lower noise of the templates, 
we can compare them with a large number of galaxies. 
Moreover, even the fits for galaxies with no spectroscopic 
redshift could be reliable in a statistical way because our 
convolutions with the SHARDS filters are interpolations 
between broad-band fluxes. Therefore, the comparison 
with the templates is another good test of the absolute 
flux calibration, and it indeed provided reassuring results 
for all filters. 

In summary, the flux calibration of the SHARDS data 
was primarily based on the comparison with HST/ACS 
grism spectra. This calibration was tested through a 
comparison with ground-based spectra and stellar popu- 
lation synthesis models fitting the broad-band photom- 
etry. Based on these three different calibrators, we also 
estimated the typical uncertainty in the zeropoints of the 
SHARDS images, typically 0.05-0.08 mag (given in Ta- 
ble [J). 

3.5. Data quality 

The SHARDS data presented in this paper were ob- 
tained in queue mode, and a maximum seeing threshold 
was set to 0.9". This image quality was imposed in order 
to reach the necessary flux limits to detect and reliably 
measure absorption indices for massive quiescent galax- 
ies up to z^2.5 (26.5-27.0 mag at the 3cr level). In Ta- 
ble [I] we show the final seeing of the SHARDS images, 
measured after mosaicking all the data for a given fil- 
ter and carrying out all reduction and calibration steps. 
Virtually all SHARDS data were taken under sub-arcsec 
conditions. 

Concerning the final depth of our images, we give in 
Table Q] the typical depths reached at the 3<r level (mea- 
sured in sextractor-based mag_auto apertures), and the 
magnitude level corresponding to the third quartile of the 
brightness distribution. The desired depths were reached 
in most of the filters, with some of them presenting ex- 
ceptional seeing and depth figures. 

4. SHARDS SCIENCE VERIFICATION: EMISSION-LINE 
GALAXIES 

The analysis of emission lines in the spectra of star- 
forming galaxies (SFGs) and objects hosting AGN is 
one of the main tools to understand galaxy evolu- 
tion. Emission lines can be used to select both SFGs 
and AGN, and then to obtain estimations of relevant 
physical parameters, being redshift, SFR, metallicity, 
or b lack hole mass some of the most interesting (see , 
e.g.. iCharlot fc Longhetti|[200lt [Kewlev fc Dopital [20021: 
iGilbank et al.ll2010t lYates et al.ll20T^ T 

Although the most straight-forward way of getting 
emission-line identifications and fluxes is through spec- 
troscopic observations, these data are hard to obtain for 
large number of objects and very time consuming. More- 
over, current spectrographs on the largest telescopes are 
only able to reach continuum magnitudes around i?i~24- 



25 mag. This results in a scarcity of spectroscopic data 
for faint objects, especially at z>1.5, and a bias of spec- 
troscopic surveys towards bright emission-line galaxies at 

A powerful alternative to spectroscopic surveys can be 
found in (ultra-)deep narrow-band imaging observations, 
which have been demonstrated to be useful to select 
emission-line galaxies (ELGs) with the faintest magni- 
tudes, and measure important quantities such as equiv- 
alent widths (EWs ) and line fluxes (see , among many, 
i Teplitz et ail 119981: IKodaira et all 120031: iPascual et al.1 
2001 1 : IQuchi et al.ll2004J: IWillis fc CourbinllgqOSHLv et al 
l2007t iTakahashi et ali l2007t IVillar et all 120081 1 2011 



Shioya et al. 2008; Nilsson et al 



Guaita et alj 2010; H aves et al 



2009: ISobralet^]l2009 ; 



l2010at ILv et all 12011). 



In addition, although imaging data cannot directly pro- 
vide robust identifications and precise observed wave- 
lengths for emission-lines, they certainly help determine 
accurate photometric redshifts, even for very high red- 
shift sources, based on the detection of both emission 
and absorption features or breaks. 

4.1. ELGs at intermediate redshift 

In this section, we analyze the sensitivity and spectral 
resolution of our ultra-deep medium-band survey to se- 
lect ELGs and measure their relevant parameters (EW, 
flux). Imaging surveys in the optical aimed at selecting 
ELGs usually employ narrow-band filters, with narrow 
meaning widths around 10 nm. The SHARDS filters are 
wider, the typical FWHM is 15 nm, but the depth, pho- 
tometric accuracy, and imaging quality of the GTC data 
(see Table [J) can compensate the lower spectral reso- 
lution (R~50), compared to more classical narrow-band 
surveys. 

Figure [7] shows an example of the selection of ELGs 
with SHARDS data for one of our 24 filters, the one cen- 
tered at 687 nm (filter F687W17). The method is similar 
to that used by narrow-band surveys such as the ones 
referenced above: the flux in a given photometric band 
at wavelength Au nc is compared with the average flux 
around that wavelength. Sources with emission-lines ly- 
ing inside the central filter would present an excess of flux 
compared to the average around it, with the latter cor- 
responding to the spectrum continuum next to the line. 
In this particular example, we would expect to identify 
a population of galaxies at z^0.80-0.84 featuring an ex- 
cess in the flux seen by the F687W17 filter due to [Oil] 
emission. Note that the redshift interval is governed by 
the width of the filter as well as by the CWL variation of 
the pass-band along the FOV. Other lines could provide 
samples at different redshifts (e.g., Ly-a at z^4.6). 

Typical narrow-band surveys use a broad-band fil 
ter to determine the co ntinuum (e.g . , IVillar et al.ll200 " 
Taka hashi et all 120071: ISobral et al.l 120091 : lHaves et aL 
2010M ISobral et al.N2012L or one or several narrow- or 
medium-band filters around a given one (|Lee et al.l2012ft . 
In our case, the contiguous spectral coverage of the opti- 
cal window allows us to obtain a continuum determina- 
tion by using the adjacent filters to a given one, or even 
several filters around the central pass-band. This mea- 
surement is very robust, since it takes into account the 
intrinsic color of each galaxy in a close region around the 
spectral region of interest, rather than providing an aver- 
age continuum value in a wide wavelength range. More- 
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over, the continuum determination is not affected by the 
emission line itself (for a typical line such as [OII]A3727), 
as would be the case for surveys using broad-band filters. 

In order to understand the scatter in Figure [3 we 
have analyzed the photometric uncertainties of our cata- 
logs. The distribution of the typical errors of our data is 
shown with dashed lines. Sources whose emission in the 
F687W17 filter is brighter than the average for the adja- 
cent filters with more than the 3-cr confidence are located 
below the continuous orange line. This is the expected 
locus for emission-line galaxies. Note that the color to de- 
tect a line with medium-band filters such as ours ranges 
between 0.1 and 0.5 mag, approximately. For filters with 
larger widths (e.g., broad-band), emission lines such as 
the ones we are able to detect would be diluted by the 
ratio between filter widths, and thus would be very hard 
to detect with broad-band data (typically having widths 
a factor of 7-8 larger than the SHARDS filters, imply- 
ing an effect on the broad-band photometry of less than 
0.1 mag for the brightest lines). 

In order to test how the SHARDS data perform when 
trying to select ELGs, we have marked in Figure [7] with 
green symbols the galaxies with reliable spectroscopic 
rcdshifts for which the [Oil] line would lie within the 
F687W17 pass-bancEE We are able to recover more than 
90% of the spectroscopically confirmed galaxies with an 
[Oil] emission-line expected within the F687W17 filter. 
A visual inspection of the spectra for the non-selected 
objects (but with z~0.8) revealed very weak or even ab- 
sent [Oil] emission-lines. Therefore, we conclude that the 
SHARDS data are very effective in isolating emission-line 
galaxies at a similar line flux level as deep spectroscopy. 

It is also interesting to note that virtually all the spec- 
troscopically confirmed ELGs in Figure [7] are brighter 
than i?~24.5. This is the spectroscopic limit of the red- 
shift surveys carried out in the GOODS-N field, and the 
typical detection threshold for the vast majority of data 
taken with state-of-thc art spectrographs in 10-meter 
class telescopes. The SHARDS observations reach at 
least 2 magnitudes fainter, and thus open the possibility 
to reliably select and study fainter and/or higher redshift 
ELGs. 

Figure |8] shows the redshift distribution of the 
sources selected as ELGs by SHARDS with the 
F687W17 filter. The filled histogram represents galax- 
ies with spectroscopic redshifts, i.e., confirmed ELGs. 
The open histogram shows the photometric rcdshifts 
(|Perez-Gonzalez et al.l I2008D for all ELG candidates. 
These redshift distributions demonstrate that the ELGs 
selected by SHARDS are preferentially located in the ap- 
propriate redshifts corresponding to the most common 
emission lines detected in intermediate and high redshift 
galaxies. The majority of the selected sample are [Oil] 
emitters at z~0.8, but other redshift peaks are also ob- 
served where one expects other relevant lines such as 
[OIII], H/3, AGN lines such as CIII] or CIV, and also 
Lya (all marked in the figure). 

Our observational strategy was conceived to detect and 
measure absorption bands. In this regard, the redshift 

21 Note that these galaxies have been highlighted only because 
of their redshift, i.e., some of these galaxies may not present [Oil] 
emission, being their spectroscopic redshifts based on some other 
spectral feature. 
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Fig. 7. — Color-magnitude diagram showing emission-line galaxy 
candidates at ~687 nm in SHARDS pointing #1. The vertical 
axis shows the color between the F687W17 band and the aver- 
age magnitude in the adjacent SHARDS bands (F670W17 and 
F704W17) measured within a circular aperture of radius r=l". 
Sources with measured spectroscopic redshifts implying that the 
[OR] line lies within the F687W17 filter are marked in green, while 
sources whose [Oil] emission is located in the adjacent filters are 
marked in blue. Solid points show galaxies with a photo-z between 
0.7<z<0.9. The dashed orange lines depict the typical photometric 
uncertainty as a function of magnitude. The locus for ELGs with 
an emission line within the F687W17 filter detected with more than 
3-<t confidence is the region below the bottom orange continuous 
line. ELGs with emission-lines within the F670W17 or F704W17 
pass-bands should be located above the top orange continuous line. 
The vertical orange line shows the 5-cr detection threshold of the 
SHARDS survey in the F687W17 band. Only the sources detected 
by IRAC for which robust photometric redshifts were estimated 
(Perez-Gonzalez et al. 2008) are depicted. 



distribution in Figure [5] also shows peaks around the re- 
gions where the Mg+Fe absorption band (at ^280 nm) 
would move through the F687W17 passband (z~1.5), 
around the Balmer or 4000 A break (z=0. 8-1.0), and 
also where the ISM dust absorption peak at 2175 A could 
produce a differential attenuation from galaxy to galaxy 
(z~2). These absorptions mimic an emission-line in a 
color-magnitude diagram such as the one depicted in Fig- 
ure [7] (see also Figure . 

Figure [9] shows eight examples of galaxies selected 
by SHARDS as emission-line sources and counting with 
spectroscopic confirmation (shown in green in Figure [7]) . 
These sources have been selected using the F687W17 
filter as central band. The two sources on the top 
row are examples of low redshift galaxies detected be- 
cause of their [OIII] AA4959, 5007 emission (forming the 
peak at z~0.3 seen in Figure [8] and marked in blue and 
green). Comparing our photometry with the available 
Keck and HST/ACS grism spectroscopy for these galax- 
ies, we show that SHARDS is able to robustly identify 
emission lines and measure EWs and fluxes. Unlike the 
noisy spectra, our data allow to study the continuum 
(compare with the stellar population synthesis models 
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Redshift 

Fig. 8. — Spectroscopic and photometric redshift distributions of 
the sources selected as ELGs with the F687W17 filter. The filled 
histogram shows the selected galaxies with spectroscopic redshifts. 
The open histogram shows all ELG c andidates using photometric 
redshifts ( Perez-Gonzalez ct al. 2008). The expected redshifts for 
emitters with some of the most typical lines (e.g., Lyo, [Oil], [OIII], 
or He*) are marked with shadowed regions. We also mark spectral 
features such as the Mg-Fe absorption band at ~280 nm, or the 
2175 A dust absorption bump (which would imprint an absorption 
in the galaxy SED). 



sho wn in blue, which were fitte d to the broad-band data; 
see lPerez-Gonzalez et ai1l2008[ ). Using all the SHARDS 
photometric data-points and fitting them to stellar pop- 
ulation models, we have estimated this continuum level 
and, combining this with the flux measurement for the 
F687W17 filter, we have obtained line fluxes and EWs. 
Note that our method to estimate the continuum is not 
based on the adjacent bands to the one containing the 
emission line alone, but on a stellar population model fit 
to all bands (except the one containing the line). When 
comparing to the values obtained from the spectroscopy, 
we typically find smaller EWs. We interpret this system- 
atic difference as an aperture effect, given that our photo- 
metric apertures are typically 1.0-1.5" in radius, proba- 
bly larger than the typical slit widths and extraction sizes 
used in spectroscopic surveys. Note that larger apertures 
would easily dilute the emission lines in a brighter contin- 
uum, thus resulting in smaller measurements of the EW. 
Flux measurements must be more robust, and indeed our 
estimations of the line fluxes are in very good agreement 
with the spectroscopic measurements (when available). 
We will present a detailed analysis of the comparison be- 
tween fluxes obtained with photometry and spectroscopy 
in a forthcoming paper (Cava et al., in preparation). 

On the second row of Figure |H1 we show two exam- 
ples of galaxies located in the ELG locus in Figure El 
but not showing emission lines. These sources were se- 
lected due to the significant flux difference between ad- 
jacent bands around the 687 nm filter. In this case, 
this large color is due to the presence of the Balmer or 
4000 A break at approximately 680 nm (as clearly shown 



by the available ground-based spectra and stellar popu- 
lation synthesis models). These Balmer break/D(4000) 
sources form the peak at z^0.7 and the redshift tail be- 
tween the z~0.8 peak, which corresponds to [Oil] emit- 
ters, and z~1.2 seen in the redshift histogram in Fig- 
ure [8] The postage stamps show that these sources are 
typically red spheroids, where strong 4000 A breaks are 
expected. Again, our SHARDS data are able to probe 
this spectral region with high accuracy and measure ab- 
sorption indices such as the D(4000). We refer the reader 
to Section [5] for a discussion on the detection of absorp- 
tion features with SHARDS data and their use for stellar 
population synthesis models. 

The third and fourth rows in Figure [9] show [Oil] 
emitters with spectroscopic confirmation (from Keck and 
HST grism data) selected in the F687W17 band ([Oil] 
around z=0.8). These examples clearly demonstrate 
that SHARDS data are able to detect faint emission- 
line galaxies, which may be even missed by low reso- 
lution spectroscopy such as that from HST/ ACS grism 
data. This is the case for the sources in the bottom row, 
whose [Oil] emission is not measurable with HST/ ACS 
data, but was confirmed with ground-based medium- 
resolution spectroscopy. We would also like to highlight 
that our medium-band spectro-photometric data are able 
to detect and robustly measure emission-lines as faint as 
^5xl0~ 17 erg s _1 cm~ 2 . 

Finally, we also point out that the postage stamps 
shown in Figure |H] demonstrate the depth and excellent 
image quality of the SHARDS images, which are able to 
detect the majority of the faint objects seen in the ACS 
data. 

In summary, our science verification test has demon- 
strated that SHARDS is able to detect virtually all the 
spectroscopically confirmed emission-line galaxies at in- 
termediate redshift. This opens the possibility to extend 
and reach higher completeness in the samples of ELGs 
in GOODS-N, and carry out a comprehensive analysis of 
star-forming galaxies and AGN down to fainter magni- 
tudes (~26.5 mag) than spectroscopic surveys. A com- 
plete analysis of the properties of the emission-line galax- 
ies detected by SHARDS will be presented in a forthcom- 
ing paper (Cava et al. 2012, in preparation). 

4.2. ELGs at very high-z 

The ultra-deep imaging data from SHARDS could, in 
principle, allow us to detect and study sources at very 
high-z (z>3). Indeed, reaching magnitudes as faint as 
26.5-27.0 mag in all bands from 500 to 941 nm is deep 
enough to detect interesting sources s uch as Lyq emit- 
ters ( LAEs) at z> 3 and up to z~6.7 ( Kashikawa et al.l 
120061 : iBouwens et all [2001 buchi et al.ll2003L l2008h . In 
Figure I10[ we show one example of a spectroscopically 
confirmed AGN at z>3 where SHARDS data reveal the 
presence of several emission lines such as CIII]A1909 or 
CIVA1548. The SED of this AGN shows that the spectro- 
photometric data from SHARDS is completely consis- 
tent and provides a similar spectral resolution as the 
HST grism data from PEARS. Our data can be used 
to obtain very robust fluxes and flux densities for both 
emission-lines and the continuum for a source which is 
as faint as the typical spectroscopic limit of ground- 
based state-of-the-art spectrographs in 10m class tele- 
scopes (see the noisy Keck spectrum for this source). 
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Fig. 9. — SEDs and postage stamps of some representative examples of the ELG candidates selected with the SHARDS F687W17 filter 
(from Figure [7}. We show galaxies which have already been spectroscopically confirmed. On the left panel for each galaxy, we depict the 
optical SED of the source. Broad-band photometric fluxes are depicted with black (ground-based) and gray (HST) filled stars. Red symbols 
show the SHARDS data, including uncertainties and filter widths. Blue lines represent stellar populati on s y nthesis mode ls fitting the broad- 
band data, which were used to obtain an estimation of the photometric redshift and stellar mass (Perez-Gonzalez et al. 2008). The cyan 
lines show ground-based spectra for each source, when available; the spectral and absolute flux calibration explained in Section I3.ll have 
been applied. For some sources, we also depict with a yellow shaded area the available HST/ ACS grism spectrum from PEARS (including 
uncertainties). The wavelength for several typical emission-lines are marked in green. We also give measurements of EWs and fluxes for 
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Fig. 10. — Examples of SEDs constructed with SHARDS data for z>3 galaxies. Symbols and lines are as in Figure [9] For each 
source, postage stamps in several bands are also given. On the left, we show a spectroscopically confirmed AGN at z~3.5 for which the 
SHARDS spectro-photometric data provide clear detections of emission-lines such as CIII and CIV. On the right, we show a dropout 
source which has been selected a s a candidate for a LAE at z~5 with a possible companion (also a LAE at the same redshift) to the NE 
(Rodrieucz-Espinosa ct al. 2012). The SHARDS data (and the very few detections in broad-band imaging observations) are fitted with a 
stellar population synthesis model that predicts an intense Ly-o emission at around 740 nm, corresponding to z=5.1. 
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Note that SHARDS data are not yet available at wave- 
lengths below ^600 nm, but when our survey is complete 
we should be able to detect and measure fluxes for Ly-a 
in this kind of object (AGN or ELGs). 

A preliminary analysis of the SHARDS data in four fil- 
ters from 687 to 738 nm revealed a dozen dropout sources 
(in one of the four filters) with SEDs consistent with z>4 
LAEs and LBGs. Figure \W\ shows one of these LAE 
candidates at z~5, including postage stamps and a SED. 
In fact, the postage stamps show two LAE candidates, 
both appearing in the F738W17 filter (one of them be- 
ing almost undetected in the ACS bands). The galaxy 
shown in the SED is a dropout in the F721W17 filter. 
Its emission is booming in the F738W17 filter, and then 
becomes fainter again (but detected) in redder filters. 
This points out to a Lyman break located within the two 
mentioned SHARDS filters, which would imply a red- 
shift around z=5. This figure demonstrates the potential 
of the SHARDS ultra-deep medium-band survey to se- 
lect magnitude 26-27 emission-line galaxies at very high 
rcdshifts. We refer the reader to Rodriguez-Espinosa et 
al. (2 012b, in preparation; see lRodriguez-Espinosa et ah! 
12012 1 for a more detailed discussion on the detection and 
properties of the LAEs detected by SHARDS. 

5. SHARDS SCIENCE VERIFICATION: ABSORPTION 
SYSTEMS 

The main goal of the SHARDS project is to analyze 
in detail the stellar populations in massive galaxies at 
high redshift, especially those that are already evolving 
passively. The observational strategy of the survey was 
devised to be able to construct rest-frame UV/optical 
SEDs for this kind of sources with enough spectral reso- 
lution to measure absorption indices correlated with im- 
portant physical parameters such as the age of the stellar 
populations. 

The main absorption index targeted by SHARDS 
for high-z passively evolving sources is the Mg index, 
Mguv- This index probes several absorption lines (e.g., 
MgIA2852, MgIIAA2796,2804, FeIIAA2600,2606) and has 
been shown to be an extremely reliable means to detect 
quiescent massive galaxies. In addition, the absorption 
lines can be used to easily distinguish the SED of a qui- 
escent massive gal axy from the featu reless spectrum of 
a dusty starburst (|Daddi et al.ll2005al ). The relative in- 
tensity of these absorptions can be measured with the 
Mguv spectral index, directly linked to the age of the 
stellar population (e.g., Bruzua l fc Chariot] Il993) . The 
index is easily and robu stly measurable eve n in low reso- 
lution spectra fR<100. IDaddi et al.|[2Q05al) . The Mg uv 
index has been successfully used in the past to obtain 
redshifts and a ges of stellar populat i ons in massive galax- 
ies at high-z dSpinrad et all Il997t iCimatti et al.1 1200 1 
McCarthy et all I2004L ISaracco et al.l 120051: IDaddi et al 
2005aHFerreras fc Yfll2004t iFerreras et al.H2012h . 

Other interesting absorption features probed by 
SHARDS (for galaxies at different redshifts) are the 
Balmer and 4000 A break, or the Ca HK lines, the 
G-band, and the Mgi, Mg 2 and Ti02 molecular bands 
(among others). All these have bee n extensively used to 



study stellar population ages (e.g., iBurstein et al 



1984; 



1999; 



IKriek et al.l[2006l l201lT iMuzzin et al.1l20Tl see also the 
review by iRenzini I2006D . The spectral resolution of the 
SHARDS dataset is adequate for measuring the Balmer 
or D(4000) breaks with high accuracy. For galaxies with 
spectroscopic redshifts or very accurate photometric red- 
shift (errors below 1%), we have very similar accuracies 
to those achieved in standard spectroscopic studies of 
nearby galaxies and synthesis models, wh ich typically us e 
bands of 10-20 nm around the feature ([Bruzual 1983). 
In the case of the molecular bands, the widths of the 
SHARDS filters are larger than the amplitude of these 
absorptions, but, given their strength, they should have 
a measurable effect on the photometry through filters 
such as ours. In any case, our filters are narrower (by 
at least a factor of x2) than those used by other op- 
tical and NIR medium-band surveys (such as MUSYC 
or NMBS). This allows us to probe these features with 
higher resolution. In addition, as demonstrated in the 
previous Section, SHARDS data are also very sensitive 
to emission lines, so they can be used to detect low level 
residual (unobscured) star formation in massive galax- 
ies at high-z and confirm their quiescent state (jointly 
with other methods, such as the study of the MIR/FIR 
emission) . 

5.1. Study of quiescent massive galaxies at high-z 

In this Section, we demonstrate the power of the 
SHARDS dataset to accurately measure absorption in- 
dices and study high-z massive galaxies. To do so, we 
focus on quiescent massive galaxies with reliable spec- 
troscopic redshifts above z=l. Rather than measuring 
and studying one specific absorption index, we com- 
bine all the photometric information for each source 
and fit the entire SED with stellar population synthe- 
sis models, i.e., we include SHARDS fluxes and also 
broad-band data up to the MIR wavelengths probed by 
Spitzer/IRAC. We, therefore, use all the (photometric) 
observations available for our sample of quiescent massive 
galaxies to carry out the most reliable analysis possible of 
their stellar population properties. We concentrate our 
study on z^l-2, a critical redshift range in which mas- 



of their mass (lArnouts et alJ 2007: 


Perez-Gonzalez et al. 


2008 


; Ilbcrt et al. 2010; Nicol et al. 


2011; Brammcr et al. 


2011 


). In this redshift interval, the D(4000) and/or 



Mguv absorption features are probed by the SHARDS 
data. 

The sources compiled for this science demonstration 
exercise have been extracted from several sources in the 
literature. The sample includes: 

• the three galaxies classified as red nuggets whose 
dynam ical masses were studied in iNewman et al.1 
(|2010f ) and are covered by our data. All the sources 
in this paper within our surveyed area were de- 
tected in all SHARDS individual filters, but some 
of the other sources from Newman et al.'s sources 
lie outside of our FOV. 



Worthevl 119941: iWorthev et al.l|1994l:jBalogh et al 
Lotz et all 120001: iFranx et al.| [20031: iKa uffman n et al 
2003bi: lForster Schreiber et all 2004; Gall azzi et al.ll200 



the four galaxies at z>l from lFerreras et all ((2009), 
who concentrated their study of early type galaxies 
on intermediate redshift sources, but had a few z>l 
galaxies in their sample. Given that the galaxies 
in this paper count with pseudo-spectroscopic red- 
shifts from HST optical grism data, we compared 
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the quoted redshifts with ou r photometric redshifts 
(|Perez-Gonzalez et alJ feOOS) and only selected the 
four galaxies with consistent estimations. 

• Finally, we also selected those galaxies with spec- 
troscopic redshifts that, according to their ob- 
served optical and NIR colors, qualify as Extremely 
Red Objects (ER O, f ollowing the defi nition in 
lElston et all 119881 and lYan et all I2000bl). and /or 
Distan t Red Galaxies (DRGs) at z>l iFranx et all 

pool . 

In order to avoid dusty starbursts in this study and 
concentrate only on massive galaxies that are evolv- 
ing passively, we only considered galaxies undetected in 
the MIPS 24 /xm data, which reach fluxes as low as 
30 yJy in GQODS-N (reducti on a nd catalogs described 
in | Perez-Gonzalez et al.1 120051 and IPerez-Gonzalez et al.1 
12008( 1. This sets an upper limit of 1-10 M yr" 1 at 
z=1.0-1.5. The final sample considered in the following 
discussion is composed of 27 galaxies at 1.01<z<1.43, 
with virtually all of them having masses above 10 10 5 M . 
Their main properties are given in Table [2] Note that we 
have only considered galaxies with a spectroscopic red- 
shift, so the sample is highly biased towards low redshift 
(close to unity, the limit lower of our sample, 75% of the 
sample is at z<1.3 and the average redshift is z=1.17) 
and not (necessarily) representative of the entire popu- 
lation of massive quiescent galaxies at z~l. 

5.2. Estimating stellar population properties from SED 
synthesis models 

Our goal in this Section is to improve the characteri- 
zation of the stellar populations of quiescent galaxies by 
using the unprecedented spectral resolution and depth 
provided by the SHARDS photometric observations. Re- 
cent results based on data with similar medium-band 
filters and (low resolution) grism spectroscopy have al- 
ready shown the power of probing key spectral features 
on the SEDs of quiescent galaxies, such as continuum 
breaks and absorption lines/bands, to place more ro- 
bust constraints on relevant properties of these galax- 
ies su ch as the age, stellar mass, or star-formati o n his- 
tory dFerreras et all [200l iWhitaker et al.l 120101 I20TT1; 
van Dokk um fe Bra mmcr 20101). 

Here we study the stellar properties of a sample of 27 
massive quiescent galaxies with spectroscopic redshifts 
above z=l by fitting their UV-to-MIR SEDs with stellar 
population synthesis (SPS) models. To calibrate the im- 
pact of using the SHARDS data in this exercise, we con- 
duct two kinds of tests. First, we show how including the 
SHARDS medium-band data allows us to decrease those 
uncertainties in the derived stellar parameters linked to 
well-known degeneracies (such as age-dust, or age-SFH). 
To do so, we analyze the connection between the pho- 
tometric uncertainties and the best-fit solutions in the 
parameter space formed by the relevant stellar proper- 
ties. Second, we study the different results obtained by 
fitting the SEDs with different sets of SPS models ex- 
tracted from several of the most common libraries found 
in the literature. 

We start the analysis of the SEDs of massive quies- 
cent galaxies at z^l by describing the set of libraries, 
the modeling assumptions, the parameter space, and the 



statistical approach to the analysis of the different solu- 
tions and degeneracies. 

The UV-to-MIR SEDs for the sample of 27 massive 
passively evolving galaxies at z>l were fitted with SPS 
models from several of the most common libraries found 
in the literature. We considered star formation histo- 
ries (SFHs) following general exponential decaying mod- 
els (r-models) using the predictions from the following 
codes: 

• t helBruzual fc Charlotl ()2003| ) models (BC03) with 
a iChabrierl (|2003ft IMF. This was chosen as our 
fiducial model (BC03-chab, hereafter). In this sec- 
tion and the following, we discuss the properties of 
our galaxies as inferred by other models in com- 
parison with those obtained with this library and 
IMF. 

• the Chariot & Bruzual (2009) models (CB09) with 
a Kroupa and Chabrier IMFs (CB09-krou and 
CB09-chab, hereafter). These models are thought 
to be an improvement over the BC03 library with 
a better treatment of some evolutionary stages of 
stars. 

• the iMarastonl (|2005l ) models (M05) with a Kroupa 
IMF (M05-krou). These models should also give 
us information about the effects on the derived pa- 
rameters of specific stellar evolutionary stages such 
as the thermally-pulsating AGB phase. 

• the PEGASE code dFioc fc Rocca-Volmerangd 
I1997L P01), assuming a lKroupal (|2001l ) initial mass 
function (IMF). These are referred to as POl-krou 
models. 

In all cases, we considered the following mass limits for 
the IMF: 0.1<A4<100 M Q . 

The fitting proce d ure i s ex plained in detail in 
IPerez-Gonzalez et all (j2003l ) and IPerez-G onzale z et all 
(|2008l) . Briefly, the photometry is compared with the 
models assuming an exponentially decreasing SFH char- 
acterized by a t parameter which runs from 1 Myr (al- 
most identical to a single stellar population, SSP) to 
100 Gyr (almost constant SFH). We run models corre- 
sponding to the different discrete values of the metal- 
licity given for each library, which run from i? /2OO to 
2.5Z e for BC03 and CB09 models, Z /2OO to 5Z for 
P01, and Z Q /50 to 2Z Q for M05 models. We considered 
an extinction parametrized with a U-band attenuation, 
A(V), with values ranging fr om to 2 mag , and assumed 
the extinction law of [Calzetti ct al. (2000). Finally, we 
probed all ages between 1 Myr and the age of the Uni- 
verse corresponding to the source redshift (up to ^6 Gyr 
for the closest galaxy) to search for the best fitting model 
minimizing a x 2 maximum-likelihood estimator. 

In order to study the uncertainties in the derived pa- 
rameters and to take into account the possible degenera- 
cies in the solutions, we run Montecarlo simulations for 
each galaxy. The method consisted in randomly varying 
the photometric data points with a Gaussian distribu- 
tion of width equal to the photometric uncertainty, and 
repeating the fit again with all possible models. We run 
the code 1000 times and then analyzed the set of solu- 
tions. This analysis identified the clusters of solutions 



18 



TABLE 2 

Stellar population synthesis results for red and dead galaxies at z~l. 
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BC03-chab 

CB09-chab 

CB09-krou 

M05-krou 

POl-krou 



158 
32 
40 

501 
50 



+ 18 
-16 
+5 
-5 
+5 
-5 
+51 
-54 
+5 



1 43+ 016 
4 00+ 40 

4 - uu -0 .43 
Q c n +0.33 
J - OU -0.35 



40 +0 ' 05 

u -^ u _0.04 
S0+ 05 

40+ 05 

u -^ u -0.05 

0.20tg;SS 
0.50t .°o 5 3 



1 000+ 000 
i.uuu_ 601 

1 000+ 000 

J--UUU-0 .000 

1 000+ 000 

u.ouu_ 000 

200+ 000 
u.zuu_ 000 



10.86i°;°S 
10.69l ;° 04 

10.78±r 4 

^+0.05 
-0.05 



11.05! 



11 ifi + 0.05 



16 +01S 

U ' 1D -0.14 
U.3»_„ 35 

u " 3 ' -0.34 

fl 99 + 0.25 

u - z -0.21 

2S + - 30 
U - ZO -0.26 



76.1 

85.6 
68.1 
91.3 
56.5 



CB09-chab 
CB09-krou 
POl-krou 
POl-krou 



584+ 94 
494j^5 

347 
227 



+ 135 
-45 
+ 31 



-42 



4 43+ 51 

*-*°-0.18 

4 72+ - 45 

'^-0.46 
, 77 +0.61 
J -' '-0.51 

■J- u '_0.40 



29+ 05 

uol -0.06 

5l+°- 05 
u.oi_ 05 



200 +0 00 ° 

200 +0 00 ° 

200+ 000 

200+ 000 
u.zuu_ 000 



11 06 +0 - 05 

11,UD -0 .05 

11 16+ - 04 

,0+0.05 

n.io_ 06 
11 10+ 006 

li,lu -0.06 



n qq+o.40 

- + 0.39 
^ -0.34 



0.367 



29+ 31 

u ' za -0.27 

29+ - 30 



14.4 
31.9 
18.2 
16.9 



,1123615.24+620944.52 



1.2674 BC03-chab 
CB09-chab 
CB09-krou 
M05-krou 
POl-krou 



1.0 



+ 0.1 
0.1 



637 
1261 



32 
2.0 



+ 0.3 
0.2 



1.02" t 



-, + 0.10 
-0.10 

1 02+ 010 

1 02+ 013 
1 - u ^-o.ii 

3.oo±S;l? 

4 50+ 45 



1.00+ 



-,+0.06 
-0.05 

20+ 05 

u - zu -0.05 

0.80t°;°S 
l.l0tg;g« 



200+ 000 

1 000+ 000 
'• uuu -i) .000 

1 000+ 000 

o.020±°;°°° 
0.005t°;°S« 



10.68l° ;° f 
in 77+0-05 

lu - ' ' -0.05 
,+0.04 
-0.06 



11. IV 



11 34+ 0( 

lllOq: -0.05 



0.68_ 
0.83+ 



+ 0.97 
0.46 



0.827 
0.891 



j+1.18 
-0.61 



88.9 
65.2 
67.2 
30.2 
40.9 



CB09-chab 


1 


Q + 0.1 

u -0.1 





.72; 


CB09-krou 




63+ 7 


1 


.01 


M05-krou 


1 


Q + 0.1 


4 


.00 


M05-krou 




31±| 


2 


.94; 


M05-krou 




18t 4 4 


2 


.98; 


POl-krou 


1 


°±°o.\ 


4 


.53; 


POl-krou 




20tl 


4 


.44; 



+ 0.08 
-0.07 
+ 0.12 
-0.10 
+ 0.50 
-0.35 
+ 0.33 
-0.34 
+ 0.62 
-0.26 
+ 0.44 
-0.51 
+ 0.48 
-0.39 



0.4CT 



-, + 0.04 
-0.05 
n O-.+0.05 
u - ol -0.04 

70 +0 - 06 

u - ' u -0.05 

o.8il° ; 06 

0- so±S:S; 

90+ 06 
u.yu_ ,o 4 

1 10+ 003 

1 - lu -0.05 



1.000 
1.000 
0.020 
0.020 
0.020 
0.020 
0.005 



0.000 
0.000 
+ 1.851 
0.000 
0.000 
0.000 
+ 0.000 
0.000 
+ 0.000 
0.000 
+ 0.000 
0.015 
+ 0.000 
0.000 



i n ^T+0.04 
10 - 67 -0.05 
a+0.05 
-0.04 



10.78^ 



I! 17+0.05 
' -0.04 
11 H + 0-05 
1 - L - 1 - L -0.03 
11 n + 0.06 

I - L - 1 - L -0.04 

ii -Hn+o- 04 

II- OU_ 04 

,1+0.05 
-0.04 



11. 34" 1 



85+ 1 ' 22 
0-83i ' 2 5 ° 

o.64±°-;s 

67+ ' 97 

u ' ul -0.46 

o.74i -°? 

,+0.98 
-0.55 



0.77! 



34.8 
32.8 
29.5 
29.2 
11.1 
37.2 
11.7 



J123627.17+620756.57 1.3950 



BC03-chab 

CB09-chab 

CB09-krou 

M05-krou 

POl-krou 



398± 4 4 

5oii 4 i 

63li 8 6 ° 4 



2.oot ; 2 2 5 
o.5ii° ^ 

2.60t° ; 29 9 

3 00+0- 34 
c ' UU -0.33 

n qq+o ° 9 
"■^^-0 .08 



20+ 05 
u - zu -0.05 

20+ 05 
00+ 05 

U.UU_Q QQ 

20+0- 05 
u - zu -0.05 

40+ 06 

u -^ u -0.05 



200+ 000 

1 000+ 000 

1,uuu -0 .000 
-, + 0.000 
-0.180 



0.200! 



020+ 000 

1 000+ 1 539 
r.uuu_ 603 



10.62+g;° 4 
10.3ll° o: ° 5 5 

7+0.05 
-0.05 



10.67! 



10 62 +0 - 05 

1U.DZ_ Q5 

10 71+0-05 

1U. ( l_o 04 



i-45li; 9 I 

2.28i 2 ; 9 70 

9 9S + 2.91 

l.68t?;i» 

n o, + 2.79 
Z - ZJ -1.63 



56.1 

52.5 
71.7 
49.1 
36.3 



BC03-chab 

CB09-chab 

CB09-chab 

CB09-krou 

M05-krou 

M05-krou 

POl-krou 

POl-krou 

POl-krou 



199± 20 , 
10l!i° 

22 
181 



1.03 



+ 3 
3 

+ 29 
-64 



162t 2 ^ 
997ti^ 

, + 118 



663 
496 
360 



74 
+ 50 
-51 
+ 47 



+ 0.13 
0.13 
_ +0.08 
U - ' J -0.07 



0.46 



0.50! 



-,+0.05 
-0.06 

no + 0.15 

1 >,o +0 - 20 

i - ou -0.15 

4.oii°o; 4 4 9 o 

2.5i; 9 ;- 

1 95+ ' 27 
i.yo_ 24 

1 46+ ' 19 

1 -^°_0.28 



+0.06 
0.14 

00+ 05 
u.uu_ 00 



0.30! 



-,+0.03 
-0.04 

14+ - 14 

U-11_o,09 

00+ 05 
u.uu_ 00 

29+ 05 

u - za -0.04 

22+ 40 

U - ZZ -0.05 

o-i2l „;°o 6 6 

0.32^ 9 2 3 



0.400 
1.000 
0.400 
1.000 
0.020 
0.020 
0.200 
0.400 
0.400 



+ 0.598 
0.206 
0.000 
-0.615 
+ 0.000 
0.396 
+ 0.000 
0.802 
+ 3.084 
0.000 
+ 0.000 
0.000 
0.199 
0.180 
+ 0.000 
0.200 
+ 0.602 
0.380 



10 58 +0 - 04 
iu.oo_ 04 

10 3"S+°- 04 
10 qo+O- 04 
10 46+ 05 

lu -™-0 .05 

10 41+ 06 
10 7-^+0-04 

1U. M_„ „ 5 

l0.75l ;- 
10.69« : 6 4 
10.69± o o: 06 



1 67+ 2 ' 25 

1DI -1.19 
n 17 +2.76 
-1.64 

i-90t?;S 

i-95ii: 6 3 4 8 

2 16+ 2 ' 85 
2.04j 
2.05j 

1 60+ 212 
i -" u -1.25 



4 +2.63 
-1.64 
- + 2.69 



39.1 
28.6 
15.0 
19.8 
35.4 
15.5 
28.6 
24.7 
10.4 



J123738. 71+621727. 86 



1.2907 BC03-chab 
CB09-chab 
CB09-krou 
M05-krou 
POl-krou 



r + 4 

-1+0.1 
-0.1 
500+ 55 



57! 
63l 
70l 
1.0+ 



2.62 
1.01 
1.01 
3.98 
3.47 



+ 0.29 
-0.29 
+ 0.12 
-0.10 
+ 0.10 
0.12 
+ 0.42 
0.38 
+ 0.34 
0.32 



0.20 
0.00 
0.10 
0.51 
0.20 



+0.05 
-0.04 
+0.04 
-0.00 
+0.04 
0.05 
+0.04 
0.05 
+ 0.05 
0.04 



200 +0 00 ° 

1 000 +0 00 ° 

1,uuu -0 .000 

1 000+ 000 
i.uuu_ 000 

020+ 000 
u - uzu -o.ooo 

200+ 000 



11.18! 



n+0.05 

-0.04 

10.83l , 5 5 
10.9ll° ;° 4 4 

, , o, +0.04 
ll.Jl_ „ 4 

, , 9 o + 0.05 



0.26 
1.20 
1.27 
0.46 
1.12 



+ 0.27 
-0.24 
+ 1.69 
-0.82 
+ 1.79 
0.86 
+ 0.48 
0.44 
+ 1.56 
0.74 



91.1 
100.0 

41.6 
100.0 

72.6 



CB09-krou 
CB09-krou 
POl-krou 
POl-krou 



135! 



=; + 20 
-16 

i-oiri 

257lf 
l.Olo.i 



2.65! 



- + 0.28 
-*-0.29 
, n o + 0.09 

9 cc + 0.43 

3.00l° o: f 3 



, + 0.04 
-0.04 

00+ 03 
u.uu_ 00 

+0.06 
0.05 
+ 0.05 
0.05 



o.io! 



0.29 
0.10 



200+ 000 

1 000+ 000 

200+ 000 
u.zuu_ lgl 

-, + 0.000 
-0.000 



0.200! 



n-i5l°„; o 5 5 

l0.88l° ; 9 4 4 

, , 91 +0.05 
11 -^- L -0.04 

11 21+0- 06 

11-^1-0.05 



. + 1.61 
J -0.79 
7+1.31 
-0.69 

1.24+1' 75 



1.15! 
0.97 4 



0.96 



0.85 
+ 1.61 
0.68 



36.6 
16.0 
15.9 
11.5 



Note. — Results of the stellar population synthesis modeling of the SEDs of massive quiescent galaxies at 1.0<z<1.4. Averages 
and ±1ct values are quoted for each parameter. First five rows for each galaxy present the results for the most statistically significant 
solution for each one of the SPS libraries presented in the text. Next, secondary solutions arc given (if available). (1) Name of the 
galaxy (including J2000 coordinates). (2) Spectroscopic rcdshift from the literature (sec text for references). (3) Stellar population 
models and IMF. (4) Exponential decay factor and uncertainty of the best fitting model (in Myr). (5) Age and uncertainty (in 
Gyr) of the best fitting model. (6) Extinction in the V-band and uncertainty (in mag) of the best fitting model. (7) Metallicity 
and uncertainty (solar units) of the best fitting model. (8) Stellar mass and uncertainty (solar units). (9) Median and 68% range 
for the goodness of the fit (x 2 value). (10) Statistical significance of this solution (in %); only solutions with a significance above 
10% are shown. This table includes all galaxies in the electronic version of the paper. 
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with a k-means method. Each cluster was assigned with 
a statistical significance given by the fraction of the 1000 
different solutions belonging to the cluster. Typically, 1 
to 4 clusters of solutions were identified for each galaxy 
above a statistical significance of 10%. In virtually all 
cases (all galaxies and models), one solution was dom- 
inant with more than 50% of the solution data points 
belonging to its cluster. Note that different clusters may 
have different typical x 2 values, i.e., the data can be bet- 
ter reproduced by some combination of parameters cor- 
responding to one cluster. The different solutions pro- 
vided information about the typical degeneracies of the 
study of stellar populations, such as the age-metallicity 
degeneracy or the one linked to r (SFH) and age. Tabled 
presents the different solutions (the statistically most sig- 
nificant and the secondary solutions) for each galaxy and 
the statistical weight of each of them. 

5.3. Evaluation of the stellar population synthesis 
modeling with SHARDS data 

In order to test how the spectro-photomctric data from 
SHARDS helps in the SPS modeling of high-z galaxies, 
we compare in Figure [TT] the results obtained when fit- 
ting only the broad-band data with those obtained when 
adding the SHARDS data for three representative galax- 
ies of our sample of quiescent massive galaxies at z>l. 
The details about the fits and results of the SPS mod- 
eling for these galaxies are presented in Figure [12] and 
Section Ol 

For the first galaxy, J123704.36+621335.03 
(z= 1.2420), the fits to the broad-band photometry 
alone provide a dominant solution characterized by an 
c-folding time t~200 Myr and an age of around 2 Gyr. 
There is another less statistically significant solution 
involving a shorter and younger burst (r~150 Myr and 
age 1.8 Gyr), and a negligible fraction of solutions would 
favor even lower e-folding times and ages. When adding 
the SHARDS data to the fits, a solution with r~150 Myr 
and t~1.5 Gyr shows the largest statistical weight. This 
solution is consistent, within uncertainties, with one of 
the solutions found when fitting the broad-band data 
alone, although a slightly younger burst is favored in the 
SHARDS fluxes (age of ~1.5 Gyr instead of ~2 Gyr). 

In the case of J123615. 24+620944.52 (z=1.2674), the 
broad-band data provide a wide range of roughly equally 
significant solutions. Each one has 20-30% statistical sig- 
nificance. They show different values of extinction (dif- 
ferences of up to 0.5 mag from one solution to another), 
metallicity and e-folding times (from an SSP to almost 
200 Myr), while ages do not vary wildly (all solutions 
predict ages slightly above 1 Gyr). The SHARDS data 
break most of this degeneracy and majorly favor a 1 Gyr 
old SSP with A(V)=1 mag. This is one of the largest 
extinctions in our sample, which could explain the rela- 
tively high flux at 8 /im, maybe revealing the presence 
of PAHs (see Figure IT2"j) . Moreover, this galaxy has a 
very red companion located 3" away to the NE (26 kpc), 
also at z sp =1.2630. This source is clearly detected at 
24 /im, which might be hiding some residual MIPS emis- 
sion for our galaxy, and provides an SFR~20 M Q yr _1 . A 
possible interaction between the two objects could have 
switched on the star formation and then stopped recently, 
thus explaining the relatively high extinction. 

In the case of J123627.17+620756.57 (z=1.3950), both 



the broad-band data alone and the whole photometric 
information including SHARDS data favor at the ap- 
proximately 60% significance level a SFH with an e- 
folding time around 400 Myr and ~2 Gyr old, with 
A(V)=0.2 mag. The SHARDS data help to discard some 
solutions found with only broad-band data which indi- 
cate very young and short bursts with high extinction 
(1 mag). 

5.4. Examples of the stellar population synthesis results 
including SHARDS data 

In Figure [HI we show the complete SEDs of the exam- 
ples presented in the last subsection. We also present the 
best fitting models for the different libraries considered 
in this paper, jointly with a close-up look into the rest- 
frame UV/optical range probed by the SHARDS data. 
We include several plots showing the results of the syn- 
thesis models and the clustering analysis of the solutions. 
Finally, we also provide postage stamps in the HST and 
SHARDS filters. 

The galaxy shown at the top of Figure [T3] lies at 
z^ 1.2420 and our currently available SHARDS data do 
cover the Mguv absorption feature. The SPS analy- 
sis for the BC03-chab fiducial models shows that there 
are at least 4 different clusters of solutions. The domi- 
nant cluster (76% statistical weight) is consistent with a 
short starburst (t^150 Myr) with solar metallicity and 
relatively low extinction, A(V)=0.4 mag. Among the 
other clusters, presenting all of them statistical signifi- 
cances below 10%, we find a longer, lower abundance, 
and older burst (t~600 Myr, t~4 Gyr, Z=Z Q /5), and 
very short bursts (almost SSPs) with an age between 1 
and 3 Gyr. These clusters share very similar extinction 
values [A(V)=0.4 mag], with shorter bursts being more 
attenuated. Note that the latter are significantly less 
statistically representative than the former (i.e., models 
do favor a 150 Myr burst). We estimate a stellar mass 
within an interval around io 10 - 8-111 M Q . Interestingly, 
although all different libraries obtain similarly good fits 
(i.e., similar \ 2 values), they provide significantly dif- 
ferent values for the relevant parameters, especially for 
the age and r values. POl-krou and M05-krou favor the 
longest, oldest, less metallic, and higher mass solution, 
while BC03 and CB09 models for different IMFs are all 
consistent with shorter and younger bursts. 

The second example in Figure [12] shows a z— 1.2674 
galaxy where we have also been able to measure the 
Mguv absorption (but only the reddest part) with the 
SHARDS data available so far, apart from other spec- 
tral features such as the bluest part of the 4000 A (or 
Balmer) break. The SPS degeneracies are considerably 
smaller than in the first example. Indeed, there is a dom- 
inant solution corresponding to a star-forming burst with 
r~l Myr (a SSP), age around 1 Gyr, A(V)~1.0 mag, 
sub-solar metallicity, and mass around 10 110 M . Al- 
though this solution clearly dominates (90% probability) 
for the fiducial model, BC03-chab, other libraries again 
achieve different results, although with very similar sta- 
tistical weight (i.e., one solution dominates for all dif- 
ferent models). Note that the different libraries provide 
very similar results in the UV/optical part of the SED, 
but differ significantly in the if-band and beyond the 
1.6 /im bump, where BC03 and CB09 models reproduce 
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Fig. 11. — Age and e-folding time results for the stellar population s ynt hesis modeling (using our fiducial library, BC03-chab) of the SEDs 
for the three massive quiescent galaxies at z>l presented in Figure 1121 For each galaxy, orange points show the solutions when fitting 
UV-to-MIR, broad-band data alone, while gray points show the results when adding the medium-band SHARDS data also. Average values 
for each cluster of solutions are marked with open symbols for the results using broad-band data alone, and filled circles for the solutions 
also using SHARDS data. Contours show 2<r uncertainties: red for fits to only broad-band, black for fits adding SHARDS fluxes. Statistical 
significances of each cluster of solutions are also given. Note that the scales are adjusted to each galaxy, but our modelling procedure has 
probed e-folding times from 1 Myr to 100 Gyr, so only a small part of the considered parameter space is shown. 



the data more closely. 

Finally, the last example in Figure Q2] shows one of the 
sources at the highest redshifts in our sample of mas- 
sive quiescent galaxies, z~1.4. In this case, two solutions 
show similar significances, with e-folding times r^200- 
400 Myr, age around 1-2 Gyr, A(V)=0.2-0.4 mag, sub- 
solar metallicity and mass 10 10 6 M Q . Similarly to the 
behavior of the different libraries in the previous exam- 
ple, in this case different models also obtain different 
estimations for the relevant parameters, especially the 
age. £/-band data probing the wavelength range around 
the rest-frame i?/-bands would most probably solve com- 
pletely the degeneracy between the models. 

5.5. Goodness of the fits 

It is interesting to discuss the properties of the sam- 
ple and the goodness of the fits for the different SPS 
libraries. Concerning the latter, a very significant frac- 
tion of the whole sample, around 95% of all galaxies, is 
better fitted by the BC03-chab models, i.e., these models 
provide better \ 2 values than any other code and IMF. 
This is the main reason for our choice of the BC03-chab 
models as fiducial in this paper. The rest of best fits cor- 
respond to the CB09-krou and POl-krou models (~2-4% 
of sources each). After the best fitting code, the second 
most preferred library and IMF is M05-krou, which gives 
the second best solution for 93% of the sources, with the 
rest of second best solutions corresponding to POl-krou. 
Note that among the two models sharing the same library 
but with different IMFs, CB09-chab and CB09-krou, the 
Kroupa IMF provides the best results, but since the best 
fitting code is BC03-chab, it seems that the goodness of 
the fit is more related to the SPS library than to the 
choice of IMF, i.e., based on the analysis of the goodness 
of the fits we cannot claim that one IMF works better 
than the other. 

To further test this result indicating a preference for 
the BC03 models for the vast majority of sources, we re- 
peated the SPS analysis but only using the broad-band 
data, i.e., without introducing the SHARDS fluxes. In- 



terestingly, the best fits were then achieved by the BC03 
models for 15% of the sample, while CB09 and M05 mod- 
els obtained the best fits for 20-30% of the galaxies, and 
P01 for just 8%. This clearly demonstrates than the 
UV/optical part (more precisely, the 200<A r „f<400 nm 
interval) of the SED is better reproduced by the BC03 
models, while the fits to the global broad-band pho- 
tometry up to the IRAC bands is similarly reproduced 
(roughly) by all codes (except P01). 

Concerning the statistical significance of the best so- 
lution, typical values for all models are always higher 
than 75%, although they could be as low as 40% a few 
cases. Three quarters of the sample have a dominant 
solution above 50% for all models and IMFs. The ac- 
tual individual values of the statistical significance of the 
different clusters of solutions are similar for the different 
codes, but for virtually all sources the most concentrated 
solution is given by the BC03-chab models (the typical 
probability for the best solution is 90%) and M05-krou 
(88%). For the fits to broad-band fluxes only, these fig- 
ures are smaller by ^20%. Figure [11] showed the compar- 
ison of SPS solutions obtained when using only broad- 
band data, and when adding the SHARDS medium-band 
photometry, for the sources plotted in Figure [T2l 

5.6. Statistical properties of the sample and systematics 

Statistics about the main properties of the stellar pop- 
ulations in our sample of massive quiescent galaxies are 
given in Tabled This statistical information is presented 
for the distinct SPS codes, and discussed in this section. 

According to our fiducial BC03-chab models, the me- 
dian stellar mass of our sample of 27 massive quiescent 
galaxies at z>l is 10 10 7 M . Typically, their stellar pop- 
ulation is 1.5 Gyr old and was formed in a burst with 
solar or slightly sub-solar metallicity with a SFH charac- 
terized by an e-folding time r~150 Myr, and currently 
presents a moderate extinction A(V)~0.5 mag. 

If we consider the average of the different stellar mass 
estimates obtained with the distinct codes and IMFs for 
a given galaxy, its typical rms is 0.15 dex, and for some 
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Fig. 12. — Three examples (one per row) of the stellar population synthesis models carried out on quiescent massive galaxies with 
spectroscopic redshifts at z>l. Each row shows the entire UV-to-MIR SED of a galaxy, jointly with the best fits for different modeling 
libraries (see text for details), and an RGB 10" X 10"postage stamp of the source built from flST data. The results for the best fitting 
model are also given for each library, including (Z /Zo,T/Gyr,t/Gyr,A(V),M./'MQ). Photometric data points include uncertainties and the 
width of the filter. To the right of the general SED, we show a close-up of the spectral region covered by SHARDS. On the top-right 
corner, we depict a postage stamp of the galaxy in the F687W17 SHARDS band. Its size is 10" X 10", and the l"radius red circle marks 
the galaxy. The bottom-right panel shows the position of the different model solutions derived with our Montecarlo method to determine 
uncertainties and degeneracies. The data points correspond to our fiducial model: BC03-chab. We depict all the 1000 solutions in different 
planes: age-T, age-extinction, age-metallicity, and age-mass. For the most statistically significant cluster, we mark the average solution 
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galaxies the rms can be as high as 0.27 dex. Most of this 
rms is due to systematic differences between the results 
obtained with different models and IMFs. According to 
our fiducial models, BC03-chab, masses for individual 
sources run from 10 10 3 M Q to 10 112 M Q . BC03-chab 
provide, on average, larger masses than CB09-krou by 
^20% when comparing values for each galaxy, with dif- 
ferences as large as a factor of 2 (in both directions). 
CB09-chab models give masses smaller than BC03-chab 
by 30% on average, and up to a factor of 3. We note 
that differences between the two CB09 models with dif- 
ferent IMFs are not just a constant offset, but are vari- 
able from galaxy to galaxy and dependent on other pa- 
rameters such as age. M05-krou and POl-krou provide 
heavier masses than BC03-chab by 15% and 40% on av- 
erage, respectively, with differences as high as a factor of 
3-4. Note that the effects of the TP-AGB phase, which 
could in principle be studied comparing classical models 
with libraries such as M05 or CB09, are complex, some- 
times producing larger and sometimes smaller masses. 

Extinctions show better agreement between differ- 
ent models. The average value for the sample is 
A(V)=0.5 mag for our fiducial models, with a typical 
uncertainty of 0.2 mag. This is also the typical average 
systematic offset from one model to another. All models 
provide systematically lower extinctions than the BC03- 
chab models. 

Our sample presents ages in the range 0.5 to 3.5 Gyr 
for our fiducial model BC03-chab. On average, BC03 
and CB09 models provide similar ages within 10%, but 
differences from model to model among these libraries 
range from 0.2 to 3 Gyr for individual galaxies. M05-krou 
and POl-krou models provide, on average, older ages by 
90% and 30%, respectively. These ages and systematic 
differences are in good agreement with the r esults found 
in the literature for z~l galaxie s (see, e.g. lHathi et"aT1 
2009t IWhitaker et alJl20lo[ 1201^ . 

The typical star formation time-scale for our sample is 
r~150 Myr for our fiducial models. For CB09 models and 
the two IMFs considered, r values are, on average, 25% 
shorter. M05 and P01 results indicate longer e-folding 
time-scales by 70% and 20%, respectively. 

Concerning metallicities, for the BC03 models the aver- 
age value is solar, and that is also the case for all the other 
codes, except M05 models, which provide an average 
value of OAZq. Note that a photometric study such as 
ours has severe difficulties in determining accurate stellar 
metallicities and breakin g degeneracies involving age, ex- 
tinction and metallicity (|0'Connelllll986t lWorthevlll99l 
iWuvts et all [2009;. but it is significant that all models 
provide consistent values near or slightly bel ow the solar 
value (c onsistently with oth er works such as | Hathi et aL 
120091 or iFerreras et al.ll2009t see also lDorman et al.ll2003 
and iKavirai et al.ll2007l ) . 

In summary, there is a significant scatter in the pre- 
dicted parameters for our sample of massive galaxies by 
different SPS codes and IMFs. It seems that this scatter 
is preferentially linked to the stellar population codes, 
rather than to the IMF. BC03 (and Chabricr IMF) pro- 
vide in virtually all cases the best solutions in terms of 
goodness of the fit and statistical significance of the best 
solution, although the x 2 differences between models are 
not statistically significant for individual galaxies. This 
preference is directly linked to the SHARDS data (prob- 



ing the rest- frame UV/optical range), given that it is 
no longer seen when fitting broad-band fluxes only. The 
stellar population in our sample of massive galaxies at 
z~l is typically (according to our fiducial models, BC03- 
chab) 1.5 Gyr old, with a formation time-scale around 
100-200 Myr, solar or slightly sub-solar metallicity, and 
moderate extinction A(V)~0.5 mag. 

5.7. Global picture of the formation of massive galaxies 

at z^l 

Our sample of massive quiescent galaxies at z>l is 
small and driven by the availability of spectroscopic 
redshifts. Consequently, its representativeness of the 
global population of such type of galaxies is highly un- 
certain. However, it is large enough to study whether 
there are significant differences in the stellar popu- 
lations of galaxies as a function of mass. This is 
the base for the downsizing scenario of galaxy for- 
mation. To analyze this topic, we have divided our 
sample in two mass bins with the same number of 
galaxies. Note that the median mass of our sam- 
ple is 10 10 - 7 M Q for the BC03-chab models, slightly 
smaller than the M*~10 10 - 9 M w value at z=1.0-l .f£f| 
(jPerez-Gonzalez et al.|[2008t iMarchesini et al.l 12009ft . so 
our (limited and not statistically representative) sample 
does probe the masses around the knee of the mass func- 
tion at redshift unity. 

Figure [13] shows the SFHs of all galaxies in our sam- 
ple for the SPS results obtained with the two best-fitting 
models, BC03 and M05. We have used a rainbow color 
code based on stellar mass: massive galaxies are plotted 
in yellow and redder colors, while lighter galaxies are de- 
picted in green or bluer colors. According to the BC03 
models, there is not a large differences in age and SFH 
between galaxies of different mass. However, there is a 
trend for the oldest galaxies to be among the most mas- 
sive. In the case of the M05, a segregation in mass is more 
evident: the lightest galaxies preferentially present young 
stellar populations, while most of the heaviest galaxies 
were formed 3 Gyr before the epoch of observation (z>2). 
In addition, some of the most massive galaxies present 
SFHs close to a SSP, also according to the M05 results. 

The ages for the lightest half-sample range from 0.6 to 
2.5 Gyr for BC03 models, and 0.2 to 3.5 Gyr for M05. 
For the most massive half of the sample, ages range from 
1.0 to 2.7 Gyr according to BC03, and 2.9 to 5.0 Gyr for 
M05 results. On average, the most massive half of our 
sample harbors stellar populations that are 1.6 Gyr old, 
while the lightest-half are 20% younger for the results 
obtained with BC03 models. For M05 models, the differ- 
ence is significantly larger: the heaviest galaxies present, 
on average, 4 Gyr old stellar populations, compared to 
2 Gyr old stars for the lightest systems. Other libraries 
provide intermediate results, i.e., a segregation in age of 
galaxies of different masses with larger differences than 
those for the BC03 models, but smaller than for the M05 
code. If we only use broad-band data in our fits, the 
results are roughly unchanged for the M05, CB09, and 
P01 libraries, but the segregation is more significant for 
the BC03 code: it provides an average age of 2.5 Gyr for 
the most massive galaxies, and 1.5 Gyr for the lightest. 

22 Value obtained after accounting for differences in IMF and 
stellar population library between different works. 
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TABLE 3 

Statistical properties of the stellar populations in z~1 massive quiescent galaxies. 
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For our fiducial models, most galaxies present e- folding 
times around 100-300 Myr and peak SFRs around 400- 
500 Mq yr _1 , but there are some examples which are 
characterized by shorter bursts and could reach even 
higher SFRs. For M05, e-folding times are larger 
(100-600 Myr) and peak SFRs are smaller. These 
SFR figures and duty c ycles are typical of ULIRGs 
and/or SMGs at z~2 (e.g.,IBlain et aij|199ftlSmail et all 
2002 IChapman et all [20051: [Perez-Gonzalez et all 120051: 
Dev et all .20081: iPope et~all ,20081) . Thev are also con- 



sistent with the SFRs predicted by cos mological simula- 
tions involving gas-rich mergers (e.g.. lDave eTalll2010t 
INaravanan et al.ll2010b llal). 

Our results are r oughly cons istent with tho s e desc ribed 
in iFerreras et all (120091) and IKavirai et all (|2012f) . In 
iFerreras et all (|2009D . thev used CB09 models to fit opti- 
cal low resolution spectroscopic data of intermediate red- 
shift massive early-type galaxies. Based on these models, 
they found typical ages of 2-3 Gyr for their last red- 
shift bin (z>0.9), slightly dependent on stellar mass. No 
major differences were found in this paper in the star 
formation timescale as a function of mass. This agrees 
well with our ages around 1-2 Gyr for z^l.2 (there is 
1.2 Gyr difference in look-back time between z=0.9 and 
z=1.2), and the very similar e-folding time s for galaxies 
with different masses. IKavirai et all poll , using CB09 
models fitting broad-band data, found similar ages for 
massive galaxies up to M~10 10 5 M Q , and a trend at 
higher masses where heavier galaxies are older. However, 
we stress that our study points out that the mass-age re- 
lationship should be taken with caution, given that it 
is highly dependent on the SPS code used in the SED 
fitting. 

6. SUMMARY AND CONCLUSIONS 

We have presented the basic characteristics of the 
Survey for High-z Absorption Red and Dead Sources 
(SHARDS), an ESO/GTC Large Program awarded 180 
hours of observing time with the OSIRIS instrument 
on the 10.4m Gran Telescopio Canarias (GTC). The 
SHARDS project was devised to be able to measure 
absorption indices such as the Mguv or D(4000) for 
galaxies at z=0.0-2.5 through imaging data, and detect 
emission-line galaxies up to z^6. For those purposes, 
SHARDS is obtaining imaging data in the GOODS-N 
field through 24 medium-band filters covering the wave- 
length range between 500 and 950 nm. The typical width 
of our filters is 15 nm for wavelengths bluer than 880 nm, 
and 25-35 nm for the 3 reddest filters. In each filter, we 
have shown that SHARDS is able to detect 26.5 mag 
sources at the 3<r level (at least). In addition, virtually 
all the images were obtained under sub-arcsec seeing con- 



of 27 massive quiescent galaxies at z=1.0-1.4 obtaincc 
given for all the sample and for the two half-sample, 
ch model). 



ditions. Currently, as of semester 2012A, SHARDS has 
obtained ~75% of its data. 

In this paper, we have presented the special reduction 
and calibration procedures used to analyze the SHARDS 
observations. Basically, the OSIRIS instrument provides 
data in a 7.5' x 7.5' region through a given filter, but each 
point in the field-of-view sees a different pass-band. The 
shape of the pass-band is constant throughout the whole 
FOV, but its central wavelength shifts to the blue as we 
move away from the optical axis, which is located to the 
left, outside the FOV. This translates to a shift of around 
10-15 nm for the CWL of the pass-band from one edge of 
the FOV to the other. Consequently, sky flat-field images 
are hard to obtain, especially for a project such as ours 
using medium-band filters which are narrow enough to be 
subject to the (variable in wavelength and time) effects 
of the sky emission. To account for this, we have devel- 
oped an OSIRIS/SHARDS pipeline that applies signifi- 
cant illumination corrections and performs an absolute 
photometric calibration in 2 dimensions for each filter. 
The flux calibration has been carried out by comparison 
with HST and ground-based (Keck, CAHA 3.5m tele- 
scope) spectra. The zeropoints have also been tested by 
comparing with synthetic magnitudes obtained from stel- 
lar population synthesis models fitting broad-band data. 
Overall, our procedure achieves an absolute photometric 
calibration uncertainty lower than 8% (typically 5%-7%). 

In this paper, we have carried our a science verifica- 
tion of the main goals of SHARDS. We have shown how 
the SHARDS data are able to detect emission lines for 
low, intermediate and high redshift sources. We have 
demonstrated that the depth and image quality of our 
survey allow us to select virtually all emission-line galax- 
ies which have been already confirmed by the deepest 
spectroscopic surveys carried out in the GOODS-N field, 
such as TKRS and DEEP3. We have shown that we are 
able to extend these spectroscopic studies of emission- 
line galaxies to faint magnitudes (fainter than the typi- 
cal spectroscopic limit, J?/~24-25 mag) and detect, for 
example, Ly-a emitters at z~5. By combining all the 
SHARDS data and fitting stellar population synthesis 
models, we have also shown that we can measure equiva- 
lent widths and fluxes for those emission lines. A compre- 
hensive study of the population of emission-line galaxies 
detected by SHARDS in GOODS-N will be presented in 
a future paper (Cava et al. 2012, in preparation). 

We have assessed the ability of the SHARDS data to 
detect, measure and model absorption features in inter- 
mediate and high-z galaxies. For that purpose, we have 
analyzed the SHARDS data for all passively evolving 
galaxies at z>l counting with a spectroscopic redshift, 
whose SED should be dominated by absorption features. 
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Fig. 13. — Star formation histories for 27 spectroscopically confirmed quiescent massive galaxies in GOODS-N. The rainbow-based color 
code differentiates galaxies according to stellar mass, with the most massive sources plotted in yellow or redder colors, and less massive 
galaxies plotted in green or bluer colors. The left panel shows the results for the BC03-chab stellar population synthesis models, and the 
right plot shows the results for M05-krou. 



We have constructed a sample of 27 quiescent massive 
galaxies with spectroscopic redshifts at 1.01<z<1.43 ex- 
tracted from several sources and covered by the SHARDS 
data (and detected in all SHARDS bands, although this 
was not a requirement of the sample selection). For this 
sample, we have analyzed the entire spectral energy dis- 
tribution from the UV to the MIR, including broad-band 
data from ground-based telescopes, HST, and Spitzer, 
as well as the SHARDS photometry. This analysis con- 
sists in fitting the SEDs with stellar population synthe- 
sis models assuming a exponentially decreasing SFH (r 
models) and an extinction law described by Calzetti's 
law. We have not imposed any a priori limitations on 
the metallicity, r value or age. In order to study the 
robustness of the derived parameters, we have used a va- 
riety of stellar population libraries, namely: BC03, M05, 
CB09, and P01. Our fitting method includes a procedure 
to probe the typical degeneracies in this kind of study, 
such as those linked to age-metallicity, or age-extinction. 
This procedure is based on a Montecarlo method which 
assigns statistical significances to the different possible 
solutions found by the synthesis models. 

Our results demonstrate that including the SHARDS 
spectro-photometric data in the analysis of the stellar 
populations in high-z passively evolving galaxies helps to 
break these degeneracies, although we do not fully suc- 
ceed in this task. When using the SHARDS data, the 
statistical significance of the best solution (that with the 
highest significance) is typically 10-20% larger than when 
fitting broad-band data alone. Typical values of the sta- 
tistical significance of the best solution when fitting all 
available photometry is 75% (i.e., one solution is favored 
at the 75% probability level). Interestingly, among the 
different stellar population synthesis libraries used in this 
paper, BC03 models (with a Chabrier IMF) provide the 
best fits (the ones with the best x 2 values) for more than 



90% of the sample when fitting SHARDS' and broad- 
band data. When only fitting broad-band data, all mod- 
els are (roughly) equally good in reproducing the data. 
We conclude that BC03 models are the most suited to 
study the rest- frame UV/optical part (up to ^400 nm) 
of the SEDs for our sample. 

Taking into account all different models, we find that 
the stellar populations in our sample of passively evolv- 
ing massive sources at z>l are typically 1.5-2.0 Gyr old, 
presenting typical r values around 100-200 Myr, solar 
or slightly sub-solar metallicity, and U-band extinctions 
around 0.5 mag. However, significant systematic differ- 
ences from one model library to another are found for 
individual sources, especially in the age and r values. 
Stellar masses are better behaved, with the largest differ- 
ences from one model to another being below 0.3dex. Us- 
ing an average of all the models, the typical masses of the 
galaxies in our sample are in the interval io 10 - 5-11 - 5 M . 

Dividing our sample of massive quiescent galaxies at 
z>l in half according to their stellar masses, we find 
that the most massive galaxies tend to be older for all 
libraries. The differences are small for BC03 results, and 
considerably larger for fits using M05, CB09 or P01 mod- 
els. Just using broad-band data also produces a clearer 
age segregation in mass. Typically, for the M05 models, 
the lightest galaxies present ages around 2 Gyr, while the 
heaviest systems are 3-4 Gyr old on average, with some 
galaxies reaching ages around 5 Gyr (i.e., they started 
forming less than 1 Gyr after the Big Bang). No signif- 
icant difference in the SFH is found between the most 
massive and lighter samples. 
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